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Abstract 

We show how to compute the probabihties of various connection topologies for 

^ ' uniformly random spanning trees on graphs embedded in surfaces. As an application, 

we show how to compute the "intensity" of the loop-erased random walk in 1? , that is, 

the probability that the walk from (0, 0) to oo passes through a given vertex or edge. 

For example, the probability that it passes through (1,0) is 5/16; this confirms a 15- 

^y • year old conjecture about the stationary sandpile density on 7? . We do the analogous 

n . , computation for the triangular lattice, honeycomb lattice and Z x M, for which the 

(-H I probabilities are 5/18, 13/36, and 1/4 — l/vr^ respectively. 
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1 Introduction 

1.1 Response matrices and groves 

Let ^ be a graph and c : ^ — t- ]R>o a positive conductance on each edge. Let A/" be a subset 
of its vertices. The triple (^,c, A/") is a resistor network. Associated to this data is the 
Dirichlet-to-Neumann matrix (also called response matrix) L, defined as follows. Given a 
function / : A^ — )■ M find its harmonic extension h on Q , that is a function on the vertices 
of Q that is harmonic on ^ \ A/" and has values / on A/". Then L{f) = —A{h)\j^ is a linear 
function of /, where A is the (positive semidefinite) graph Laplacian. In electrical terms, 
L{f) gives the current fiow into the nodes A^ when they are held at / volts. While it is not 
obvious from this definition, L is a symmetric negative semidefinite matrix. 

Planar resistor networks where A/" is a subset of the outer face were called circular 
planar networks in |CIM98j . and were studied in |CdV94l ICdVGV96l ICIM98J . These 
authors classified which matrices occur as response matrices: they are precisely the matrices 
whose "non-interlaced" minors are nonnegative. (A non-interlaced minor is one in which 
there are no 4 indices a < b < c < d for which a and c are rows and b and d are columns 
or vice versa.) Furthermore they showed how to construct a graph having a given matrix L. 

On a circular planar graph a grove is a spanning forest (set of edges with no cycles) in 
which every component contains at least one vertex in A/". In |KWllaj we studied the natural 
probability measure on groves (where each grove occurs with probability proportional to the 
product of its edge weights), showing how to compute the probability that a random grove 
has a given connection topology in terms of the entries in L . 

1.2 Graphs on surfaces 

We study here the same problem for a graph Q embedded on a surface S. Here the usual 
notion of response matrix is not rich enough to extract information about the underlying 
topological structure of a grove. Given a resistor network on a surface S, a natural general- 
ization of the response matrix is a matrix-valued function =Sf on the representation variety 
Hom(7ri(S), iif) of flat if -connections on Q; here H = C* or SL2(C). We show here how 
=Sf can be used to compute connection probabilities of (certain types of) groves on Q . 

The question of characterizing which matrices Jif occur as a function of the topology 
of S remains open. See Lam and Pylyavskyy [LPlla] for related work in the case when the 
surface S is an annulus. 

We give special attention to the case where the surface S is an annulus; this is the 
easiest case beyond the planar one (but already quite involved) and also has applications to 
the study of spanning trees on planar graphs. 



1.3 Applications to planar graphs 

Using these techniques one can in principle compute the probabihty that the path of the 
uniform spanning tree from a to 6 in a planar graph passes through a given set of edges or 



vertices (as in Figure!). We carry out this computation for Z for a single edge or vertex 



see 



Figure 2); our method shows that the answer is in Q(vr); we conjecture it to in fact be 



in Q[-] . For the triangular and hexagonal lattices, we conjecture these probabilities to be 



in Oi^l (see 



Figure 13 and Figure 15). 



For example, we show that the probability that the loop-erased walk in 1? from (0, 0) to 
oo contains the point (1,0) is 5/16. (See [Figure 1[ ) This number was predicted by Levine 



and Peres |LPllb] and by Poghosyan and Priezzhev [PPllj . by relating this probability to 
the average stationary density of the abelian sandpile model. 

The connection between the spanning trees and the abelian sandpile model was discovered 
by Majumdar and Dhar [MD92J . and Priezzhev |Pri94j used this connection to compute 
the height distribution of the abelian sandpile model, in terms of two integrals that could 
not be evaluated in closed form. Grassberger evaluated these integrals numerically, and 
conjectured that the stationary density of the sandpile on 1? is 17/8. Later Jeng, Piroux, 
and Ruelle [JPR06J showed how to express one of these two integrals in terms of the other, 
and determined the sandpile height distribution in closed form, under the assumption that 
the remaining integral, which numerically is 0.5 ± 10^^^, is exactly 1/2. Our derivation of 
this probability that LERW passes through (1,0) confirms these conjectures (although our 
methods are different), and shows that this aforementioned integral is exactly 1/2. 

Similar computations compute this probability on the triangular grid and honeycomb 
grid; the probabilities are 5/18 and 13/36, respectively. We do not know why these are 
rational numbers. On Z x M (with exponential horizontal jump rates) the probability is 
l/4-l/7r2. 

We can also compute these probabilities for other vertices. For example, in the square 
lattice, the probability that the LERW from (0,0) to oo passes through (1,1) is 1/4 — 
l/(47r) + 1/(27?^), and the probability of that it passes through (2,0) is 1/8 + l/(47r) + 
l/(47r2) - 3/(27r3) + l/(27r^) (see [Figure 2]). 



While we were writing up our results, Poghosyan, Priezzhev, and Ruelle independently 
found another proof that the probability of visiting (1, 0) is 5/16 [PPRllj . They also asked 
about the probability about visiting other points, and remarked that the probability that the 
LERW visits (1, 1) is numerically close to 2/9. (This differs from the true value by about 

10-3.) 

There are some intersting coincidences in the (undirected) edge intensities of loop-erased 
random walk. For each of the square, triangular, and hexagonal lattices, there are several 
groups of edges which are unrelated by any symmetry of the lattice for which the undirected 
edge intensities are identical (see Figures [HI [T31 and[T5]). We do not have an explanation 
for this phenomenon. 



_F 



F 



Figure 1: A portion of the uniform spanning tree on 1? (actually a portion of the UST on 
a very large grid), with the path from (0,0) to oo shown in bold. The uniform spanning 
tree on H? can be constructed as a weak limit of uniform spanning trees on large boxes. 
It is known that the limiting measure exists, is unique, and is supported on trees of 1? 
|Pem91] . (This is in contrast to e.g., Z^, where the limiting measure is supported on forests 
with infinitely many trees |Pem91] .) Almost surely, within the uniform spanning tree of Z^, 
each vertex has a unique infinite path starting from it |Pem91] (see also |BLPSOl] ILPllc] for 
modern proofs). The path to infinity (or in the case of a finite graph, the path to a boundary 
vertex) is a loop-erased random walk (LERW) |Pem91] (see also |Wil96] ). In this example, 
the LERW from vertex (0,0) to oo passed through the vertex (1,0). The probability that 
this happens is the LERW intensity at (1,0), which is 5/16. 
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Figure 2: Intensity of loop-erased random walk on Z^. The origin is at the lower- left, and 
directed edge-intensities as well as vertex-intensities of the LERW are shown. (See also 
Figure iTj ) 
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2 Bundles and connections 

Let ^ be a graph. Given a fixed vector space V , a ^-bundle, or simply a vector bundle 
i? on ^ is the choice of a vector space V^ isomorphic to V for every vertex v oi Q . A vector 
bundle can be identified with the vector space Vg = (BvVy = V'^' . A section of a vector 
bundle B is an element of Vg . 

li H is a subgroup of Aut(y), an iJ-connection $ is the choice for each directed 
edge e = {v, w) of Q of an isomorphism 0„ ,^ G H between the corresponding vector spaces 
<t>v,w '■ Vv -^ Vw, with the property that 0„^^ = (p~\- This isomorphism is called the parallel 
transport of vectors in V^ to vectors in V^,. Given an oriented cycle 7 in ^ starting at v , 
the monodromy of the connection is the element of Aut(K,) which is the product of these 
isomorphisms around 7. Monodromies starting at different vertices on 7 are conjugate. 

Two connections $ = {4>e} and $' = {0e} are gauge equivalent if there are maps 
ipv '■ Vv —> K such that (f)v,w ° i^v = "ipw ° <P'v w ^^ ^^^ vertices v and w of Q . A connection is 
trivializable if it is gauge equivalent to a trivial bundle (in which all (pv^w are the identity). 

It is useful to extend the notion of bundle and connection to the edges as well: define for 
each edge e of ^ a copy Vg of K, and define maps (f)v,e : K — > V'e whenever v is an endpoint 
of e, with the property that (f)e,w ° 4>v,e = 4>v,w whenever edge e joins vertices v and w. 

A line bundle is a V -bundle where V = C, the 1 -dimensional complex vector space. 
In this case if we choose a basis for each C then the parallel transport is just multiplication 
by an element of C* = C \ {0}. Furthermore, the monodromy of a cycle is in C* and does 
not depend on the vertex of the cycle. 

In this paper we will take K = C^ or C^, and use H = C*- or H = SL2(C) -connections. 
In some applications it is more appropriate to use U(l) C C* and SU(2) C SL2(C) as 
structure groups. 

2.1 Laplacian A 

Let c: E ^ ]R>o be a conductance associated to each edge. We then define A : Vg — )■ Vg 
acting on sections by the formula 

A/(^)= Yl c.,«,(/(t;)-0»,./H). (1) 

w:(v,w)£E 

A section is said to be harmonic if it is in the kernel of A. 

We define an operator d from sections of the bundle over vertices to sections over the 
edges, for an oriented edge e = {v,w), by 

df{e) = C,,,aj{<Pv,ef{v) - (t)n,,ef{w)), 

and its "adjoint" 

d*u{v) = y^0e,t,^(e). 
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Then the Laplacian can be written A = d*d [KenlOj . 

A cycle-rooted spanning forest (CRSF) in a graph is a set of edges each of whose 
components contains a unique cycle, that is, has as many vertices as edges. A component of 
a CRSF is called a cycle-rooted tree (CRT). 

Theorem 2.1 ( |For93t [KenlOj ). For a C* -connection, 

detA= y: n ^(^)n (^-^^-^ 

CRSFs edges e cycles ^ 

where the sum is over cycle-rooted spanning forests, where the first product is over all edges 
of the CRSF and the second is over cycles of the CRSF, and w is the monodromy of the 
cycle defined above. 

For a U(l) -connection, A is Hermitian and positive semidefinite [KenlOj . The mon- 
odromy of a loop is in U(l) and so 2 — w — 1/w > and thus we can define a probability 
measure on CRSFs where a CRSF has a probability proportional to Yle ^i^) 11(2 ~ "^ ~ — ) • 

A similar result holds for a SL2(C) -connection. Now A is a quaternion-Hermitian matrix, 

that is, a matrix with entries in GL2(C) which satisfies A^j = A*j, where 



c d 



d -h 
—c a 



. Its q-determinant counts CRSFs: 
Theorem 2.2 ( |KenlO] ). For an SL2{C) -connection, 

qdetA= J2 n ^(^) n(2-Trw), 

CRSFs edges e cycles 

where the sum is over cycle-rooted spanning forests, where the first product is over all edges 
of the CRSF and the second is over cycles of the CRSF, and w is the monodromy of the 
cycle. 

For information on q-determinants, see |Dys70| ; for the purposes of this paper one can 

define qdet M = V det M , where if M is an A^ x A^ matrix with entries in GL2(C) then M 
is the 2N x 2N matrix with C entries obtained by replacing each entry of M by its 2 x 2 
block of complex numbers. In the cases of primary interest M is a quaternion-Hermitian, 
or "self-dual", matrix; for self-dual matrices qdet is a polynomial in the matrix entries. 
Matrices with GL2(C) entries enjoy many of the properties of usual matrices: for example, 
multiplication and addition work the same way. The inverse of a self-dual matrix is well- 
defined and is both a left- and right-inverse, see e.g., |Dys70] . 

For an SU(2) -connection, 2 — Tr w > and this allows us to define a measure on CRSFs 
as above. 



2.2 Dirichlet boundary conditions 

If -B C ^ is a set of vertices, the Laplacian with Dirichlet boundary conditions at B is defined 
on sections over Q \ B hj the same formula ([I]) above with v E Q \ B and the sum over 
all Q . In other words A is just a submatrix of the usual Laplacian on Q . Its determinant 
also has an interpretation. A CRSF on a graph with boundary i? is a set of edges such that 
each component is either a CRT not containing any vertex of i? or a tree containing a single 
vertex oi B . In this setting Theorems 12.11 and 12.21 have the same statements (where tree 
components do not have any monodromy term). See [KenlOj . 

2.3 Green's function ^ 

The Green's function for the standard Laplacian (with Dirichlet boundary conditions) is the 
inverse of the Laplacian. It has the probabilistic interpretation that Gpg is (^^ Cq,r)~^ times 
the expected number of visits to g of a simple random walk started at p (and stopped at the 
boundary); equivalently, it is (^^ Cq_r)~^ times the sum over all paths from p to q which do 
not hit the boundary, of the probability of the path. 

In the case of a graph with connection, the Green's function W is again the inverse of 
the Laplacian, and has a similar probabilistic interpretation: 

Lemma 2.3. $^p_g is C^j.Cq^r)~^ times the sum over all from, p to q of the product of the 
parallel transports along the path (from q to p) times the path probability (from p to q). 
This will be matrix-valued for an SL2{'C) -bundle. 

Proof. Using the above definition of ^ as a sum over paths, 

/ ^ '^p,r^r,q '^P,q / j ^q,i' / ^ '^p,rCr,q(Pq,r 
r rr^q rr^q 

and since any nontrivial path to q must have last step from a neighbor of q , this equals zero 
unless p = q and the path has length 0, in which case the second sum is zero and the first 

term is (^^ Cq^r)~^ Y^r (^r,q = 1 • □ 

2.4 Response matrix ^ 

Let A/" be a nonempty set of nodes of Q , and n = \f/\ . For each node v pick a preferred 
basis for V^ , the vector space over v . 

We define an nxn matrix ^ = ^$ (with entries in H), the response matrix, or Dirichlet- 
to-Neumann matrix, from this data: =Sf : V — )■ V is (minus) the Schur reduction of the 
Laplacian A to N' . That is, .if is defined as follows. Order the vertices so that A/" comes 
first. In this ordering the Laplacian is 

'A B 
B* C ■ 
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A 



Then =Sf = —A + BC^^B* . Note that C is the Laplacian with Dirichlet boundary conditions 
at A^, and since Q is connected, C is (for generic connections) nonsingular. 

From the viewpoint of harmonic functions, =Sf is the Dirichlet-to- Neumann matrix: given 
/ G y , find the unique section h with boundary values / at the nodes and harmonic at 
the interior (non-node) vertices. Then =Sf/ is —Ah evaluated at the nodes. To see this, let 

M. 



hi be h at the interior vertices, that is, h 



If 



A/i 



A B 
B* C 



then B*f + Chi = 0, i.e., hi = -C~^B*f, and then c = Af - BC-^B*f = -if/. 

The response matrix ^ has entries which are functions of the parallel transports. See 
ITheorem 4.21 below for an explicit probabilistic interpretation of the entries oi ^ . In order 
to define =Sf as a matrix one must choose a basis in Vy for each node v. Base changes in 
the Vy then act on Sf by conjugation by diagonal matrices with entries in H . 

Lemma 2.4. When the Laplacian A is nonsingular, the response matrix S£ is given by 

^ = -{A-%)-\ 



A B 
B* C 
is invertible since it is the 



Proof. Write A 



where A is the submatrix indexed by the nodes. Submatrix C 
/aplacian with Dirichlet boundary conditions. We have 



A 



A-BC-^B* BC-^ 
/ 



/ 
B* C 



and using ^ = -A + BC'^B* 

A-i = 



/ 



/ 



* 
-1 



D 



If ^ is the Green's function with boundary at node n, then ^ = A~^ , where the Dirichlet 
Laplacian A is obtained from A simply by removing row and column n. Since =Sf has the 
same response as A on the set A/", the response matrix of A is just S^ with row and column 
n removed. Thus \^i,o^^i"''^i = —A^^W\{n}, or equivalently. 



[^^. 



■\j=l,...,n—l 
jJi=l,...,n-l 



'«.iJi=l,...,n-l 



-1 



(2) 



By perturbing A , we see that ([2]) holds even if the Laplacian A is singular, so long as the 
Dirichlet Laplacian is nonsingular. 
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3 Graphs on surfaces 

Let S be an oriented surface, possibly with boundary, and Q a graph embedded on S in 
such a way that complementary components (the connected components of the surface after 
it is cut along the edges of Q) are contractible or peripheral annuli (that is, an annular 
neighborhood of a boundary component). We call the pair (^, S) a surface graph. 







— 






— , 

^@ — 

i / i „ 




/ „ „ 






-^, ,, ,, ^K^, u . 




— 




Figure 3: On the left is a graph Q , which is the square grid restricted to a finite box with 
wired boundary conditions: the outer boundary is one vertex (and the bottom edge is a 
self-loop). This graph will be relevant to our later loop-erased random walk computations. 
The middle diagram shows Q embedded in a surface S, which in this case is an annulus 
whose inner boundary is one of the squares of the grid. There is a "zipper" (edges crossing a 
dual path) connecting the inner boundary to the outer boundary of the annulus, and edges 
crossing the zipper have parallel transport z from their left endpoint to their right endpoint. 
We have labeled four of the vertices on the boundary of E , which we call nodes. On the 
right is a schematic diagram of the surface graph which shows the surface E , the zipper, and 
the nodes of Q , but not the internal vertices and edges of Q . 



3.1 Nodes and interior vertices 

For each boundary component of the surface E there is a "peripheral" cycle on Q , consisting 
of those vertices and edges bounding the same complementary component as the boundary 
component. We select from this cycle a (possibly empty) set of vertices. The union of these 
special vertices over all boundary components will be the nodes M \ the non-node vertices 
are interior vertices (even though these may be on the boundary of E). 

Planar maps, in which E is topological disk, are examples of surface graphs: these are 
called circular planar graphs in |CIM98] . In this case the nodes are a subset of the vertices 
on the outer face. 
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3.2 Flat bundles 

Given a surface graph [Q, S) , a vector bundle on Q with connection $ is flat if it has trivial 
monodroniy around any loop which is contractible on S. In this case, the monodromy 
around a loop only depends on the homotopy class of the (pointed) loop in 7ri(S), and so 
the monodromy determines a representation of vri(S) into Axit{Vp) , where p is the base 
point. This representation depends on the base point p for tti; choosing a different base 
point will conjugate the representation. 

Conversely, let p G Hom(7ri(E), Aut(V^)) be a representation of 7ri(S) into Aut(V); there 
is a unique (up to gauge equivalence) flat bundle with monodromy p. It is easy to construct: 
for example start with a trivial bundle on a spanning tree of Q; for each additional edge 
the parallel transport along it is determined by the topological type of the resulting cycle 
created. 

In the case of a line bundle, Aut(C) = C* is abelian and the monodromy of a loop is 
well defined without regard to base point. Moreover in this case p{'~f) only depends on the 
homology class of the loop 7, since any map from 7ri(S) to an abelian group factors through 
i^i(S). 

4 The response matrix and probabilities 

4.1 Circular planar graphs 

In the case of a planar graph, there is no monodromy and =Sf = L is a matrix of real numbers. 
This case was analyzed by |CdV94] , see also |CdVGV96| ICIM98] . Colin de Verdiere showed 
that response matrices L of planar graphs are characterized by having nonnegative "non- 
interlaced" minors. Given two disjoint subsets of nodes R and S , we say that R and S 
are non-interlaced if R and 5* are contained in disjoint intervals in the circular order on the 
nodes. When \R\ = \S\, the corresponding minor is det(L'^) > (the determinant of the 
sub matrix whose rows are indexed by R and columns by 5). 

In [KWllat Proposition 2.8] (see also jCIMQSt Lemma 4.1]), there is an interpretation 
of the entries of L in terms of groves. A grove is a spanning forest with the property that 
every component contains at least one node. The weight of a grove is the product of the 
conductances of its edges. 

Theorem 4.1 ( |CIM98[ IKWlla] ). When i ^ j , Lij is a ratio of two terms: the numerator 
is the weighted sum of groves in which nodes i,j are in the same component and every other 
node is in a component by itself; the denominator is the weighted sum of groves in which 
every node is in its own component. 

The diagonal terms have a slightly different interpretation: Li^i has the same denominator 
as the others; hut its numerator is minus the weighted sum of components in which nodes 
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different from i are in separate components, and i is in one of these components. 

More generally, for disjoint non-interlaced subsets R, S C Af with \R\ = \S\, det(L;^) is 
a ratio of two terms: the denominator is the weighted sum of groves in which every node is 
in its own component, and the numerator is the weighted sum of groves in which the nodes 
in R are connected pairwise with nodes in S , and other nodes are in their own component. 

In particular the last statement proves that the non-interlaced minors of L are nonneg- 
ative. 

Groves can be grouped into subsets according to the way they partition the nodes (that 
is, the way the nodes are connected in a grove). For example, a grove of type 1, 2 | 3, 4, 5 | 6 
is one in which nodes 1 and 2 are in a tree, nodes 3,4,5 are in a second tree, and node 6 
is in its own tree. For a partition a of the nodes, we let Z[a] denote the weighted sum of 
groves of type a. For circular planar graphs with n nodes on the boundary, we previously 
showed |KWllaj how to compute the ratio Z[a]/Z[l\2\- ■ -In] for any planar partition a of 
{1, 2, . . . , n} . This ratio is an integer-coefficient polynomial in the Ljj |KWllaj . 

4.2 ^ matrix entries 

Like in [Theorem 4.11 in the case of a flat bundle on a surface graph there is a combinatorial 
interpretation of the entries of =Sf . A collection of edges of a surface graph (^, S) is a 
cycle-rooted grove (CRG) if each component is either a CRT (a component containing 
one cycle) not containing a node, or a tree containing at least one node. Moreover for each 
CRT component, the cycle must be topologically nontrivial. A CRG is distinguished from a 
CRSF by the fact that in a CRG the tree components may contain several nodes, while in 
a CRSF the tree components contain a unique node. 

A CRG has a weight which is the product of its edge conductances times the product 
over its cycles oi 2 — w — 1/w (for a line bundle) or 2 — Tt{w) (for a SL2(C) -bundle), where 
w is the monodromy around the cycle. For a partition a of the nodes, we deflne 

^[cr] := weighted sum of cycle-rooted groves of type a 
3f := weighted sum of cycle-rooted groves in which all nodes are connected 

For example, the weighted sum of CRSFs is j!°[l|2| ■ • ■ |^]. Suppose the partition a is a 
partial pairing, i.e., a consists of doubleton and singleton parts, say cr = ri, si | ■ ■ ■ | r^, s^ | 
ti \ ■ ■ ■ \te. We can deflne 

k 

^[rl\ ■ ■ ■ Irfcl^il ■ ■ ■ l^d ■= -^k] ^ TT parallel transport to r^ from s, 

1=1 

for line bundles (so that the structure group is commutative and the above product makes 
sense), and for vector bundles when a has only one doubleton part. 
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Theorem 4.2. If i ^ j , then 

^ [^ I (nodes other than i and j in singleton parts)\ 



^ 

=^7.. 



J^[l|2| 



\n\ 



A B 
B* C 



be the Laplacian of Q , with 



Proof. Let us first do the hne bundle case. Let A 

A indexed by the nodes M . 

We make a new graph Q by adding an edge Cij to Q which connects i and j and has 
parallel transport z when directed from i to j . Let A be the line bundle Laplacian on the 
new graph Q , with Dirichlet boundary conditions at the nodes except nodes i and j , that 



is, A 



a 
b* 



b 
C 



where a 



Ai,i + 1 



^3 + 1 



and b is the zth and jth column of B . 



By [Theorem 2.11 (and its extension discussed in section 12.21) . — [2](det A) is a sum of 
CRSFs with each node except z,j in its own tree component, and nodes z,j in a cycle 
containing edge ejj , and the weight includes the parallel transport of the path in Q from j 
to i. (Here [z°']f{z) refers to the coefficient of 2° in f{z).) We can write 



detA 

-[2]det A 



-bC-^b* bc-^ 

/ 

det[a-6C-ife1detC 




C 



(3) 



[^°][a-6C-^6*]i,2detC 



However 



-1^*- 



1,2 



[A - BC-^B*]i^j = -^i 



[z^][a-bC- 
Finally, det C is the sum of CRSFs. 

The proof in the SL2(C) -bundle case is similar. Let A 



A 
B* 



B 
C 



be the SL2 



bundle Laplacian of Q . Add an edge e, .■ to Q from node i to node j with parallel transport 



M e SL2(C). As above let A 



where a 



Ai.i + / 



Aj^i - M 



A; . - M-l' 



A, 



and b is the 



a b 
b* C 

2th and jth column of B . Now det A gives a weighted sum of CRSFs with each node except 
2, j in its own tree component, where the weight is the product of the monodromies along 
the cycles. If the CRSF contains a cycle that uses edge Ci^j , then the monodromy of this 
cycle will depend on M, and otherwise, the weight of the CRSF does not depend on M . 
We can write 

qdet A' = Co + 5^ C^(2 - Tr(ir^M)), (4) 

where the sum is over configurations u with a cycle 7 containing edge Cj j- , K^ is the parallel 
transport to i from j in the cycle 7, and Co,C^ and K^ do not depend on M . 
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We have ([3]) in this case as well, where C does not depend on M . Letting D = hC ^b* . 
which does not depend on M , we can write 



a - bC-'b 



-It* 






^i,j 



D, 



*j 



M 



-1' 



This is a 2 X 2 matrix with entries in GL2(C), whose qdet equals 

qdet(a - bC'^b*) = TtHAj - Dij)M) + Ci 
where Ci does not depend on M . Comparing with (jlj) we see that, since M was arbitrary, 
- J2 C'c.^.. = Aj - A, = [A- BC-'B%, det C = -^i,. n 



j«j 



Principal minors of ^ also have probabilistic interpretations: 

Theorem 4.3. Let S C A/", then the S x S minor det^Sf^* is the ratio of two terms: the 
denominator is the weighted sum of CRSFs; the numerator is (—1)''^' times the weighted 
sum of CRGs of Qs, the graph Q in which all nodes in S are considered interior, and in 
which the remaining nodes are in separate components. If Q = {qi, . . . ,qi} = Af \S , then 



det^^c? 



■1) 



\s\ ^[lM---\qi] 



^[1|2| 



\n\ 



Proof. Order the vertices of G by first N'\S , then S , then the internal nodes. In this order 
we have 



Ai A2 


B{ 


A\ A3 


B2 


Bl Bl 


C 



Then 



^ = - 



Ai A2 

A\ A3 



+ 



52 



C-^ {Bl B*\ 



and ^/ = det(— y43 + B^C ^B1^. The proof follows from the identity 



■^3 B2 

B* C 



A, 



D /^ — 1 D* D /^ — 1 
JD2^ -D2 D2Ly 

/ 



Bl C 



upon taking determinants: the left-hand side determinant is the weighted sum of CRGs of 
Qs^ the right-hand side determinant is (—1)'"^' times the product of ^^ and detC which 

X yA \ X 0" 

counts CRSFs. In the SL2(C)-case we used the fact that qdet n t ~ ^f^^^ v t ~ 
qdetX. D 
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For line bundles there is an interpretation of more general minors: 

Theorem 4.4. Let Q, R, S, T be disjoint sequences of M for which M = QU RU S UT and 
\R\ = \S\. Then det.SfJ'u^ is the ratio of two terms: the denominator is the weighted sum 
of CRSFs; the numerator is a signed weighted sum of CRGs of Qt , the graph Q in which all 
nodes in T are considered interior, the nodes in Q are in singleton parts, and in which nodes 
in R are paired with nodes in S , with the sign being the sign of the pairing permutation, 
times the parallel transports from S to R: 

^rsp(i)| |*p('=)|„ I |„ 1 
detJ^«,T- 2^ I iJ ^[l|2|---|n] 

permutations p l i i i j 

Proof. The previous proof shows how to reduce to the case T is empty. Furthermore 
det J^^ det C = det A^^j where I is the set of internal nodes. So we need to evaluate 
det A^^j. The proof now follows the proof of [Theorem 2. II which is found in |KenlOt proof 
of Theorem 1]. Write A = d*d where d is the operator from sections over the vertices to 
sections over the edges. Then A'^^j = d*gyjjdFtui where dx is the restriction of d to sections 
over X . By the Cauchy-Binet theorem, 



det d*suidRui = J^ det{dl^j)* det d^yjj, 



where the sum is over collections of edges Y of cardinality I^U/I . The nonzero terms in the 
sum are collections of edges in which each component is a CRT if we glue Tj to Sj for each i. 
Equivalently, each component is either a CRT or a tree containing a unique r E R and 
s G 5*. The weight of a component with a cycle is 2 — w — 1/w where w is the monodromy 
of the cycle; the weight of a path is the parallel transport to the r from the s. It remains 
to compute the signature of each configuration. 

This signature is the same as the signature in the case of a trivial bundle, which is 
determined by JCIM98J to be the signature of the permutation from R to S determined by 
the pairing. D 

4.3 Cycle-rooted grove probabilities 

Theorem 4.5. The probability of any topological type of ORG involving only two-node con- 
nections and loops is a function of S£ and the weighted sum of CRSFs. 

Proof. By a result of Kenyon |KenlO] (based on a theorem of Fock and Goncharov |FG06j ). 
on a graph embedded on a surface with no nodes one can compute the probability of any 
topological type of CRSF (that is, the probability that a CRSF has a given set of homotopi- 
cally nontrivial cycles up to isotopy) from the determinant of the Laplacian considered as a 
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function on the space of flat SL2(C) -connections. Indeed, as X runs over all possible "finite 
laminations", that is, isotopy classes of collections of finite, pairwise disjoint, topologically 
nontrivial simple loops on the surface, the products ncyciesmx(2 ~ Triy) form a basis for 
a vector space (the vector space of regular functions on the representation variety) and the 
Laplacian determinant is an element of this vector space. In other words. [Theorem 2.2l above 
shows that det A = J2x ^^ Ilcycies in x(2^Tr w) , where X runs over finite laminations; such 
an expression determines each coefficient Cx uniquely. 

To compute the probability that a random CRG y on S has a fixed topology of node 
connections, add edges Cij to Q connecting endpoints of all two-node connections i ^ j 
of Y ; the resulting graph Q' can be embedded on a surface S' containing E , and the union 
of Y and the new edges is a CRSF on Q' . (We obtain E' from S by gluing a single strip 
running from i to j for each Cj j- ; in this way cycles containing different Cij s are in different 
homotopy classes.) 

Any CRG on Q with the same node connection type as Y can be completed to a CRSF 
on S' by adding the edges Cij. Conversely (since each added edge e^j- is in a different 
homotopy class on S') each CRSF of this topological type comes from a CRG on Q with 
the same connection type as Y . 

The flat connection on Q can be extended to a flat connection on Q' by taking generic 
parallel transports (pij along the Cj j . 

It remains to show that the Laplacian determinant of Q' is a function of =Sf , the new 
parallel transports $ = {(pij} , and detC. However A^ = A + 5$ where S<s> is supported 



on the nodes; using A 



B* C 



we have 



det A^' = det 



A + 5$ B 
B* C 



det 



A + S^-BC-'B* BC 




det(-^ + 5$)detC. 



Y-i' 



/ 
B* C 



U 



5 Basic surface graphs 



The simplest non-circular-planar case is the annulus. Since vri(S) is abelian in this case 
it usually suffices to consider a line bundle rather than a two-dimensional bundle. The =Sf 
matrix then depends on a single variable z G C* which is the monodromy of a flat connection. 
For simplicity we choose a connection which is the identity on all edges except for the edges 
crossing a "zipper" , that is, a dual path connecting the boundaries; these edges have parallel 
transport z. 

Suppose (^, S) is a surface graph on an annulus with rii nodes on one boundary compo- 
nent and n2 on the other. Then =Sf is an (rii -|- 712) -dimensional matrix with entries which 
are rational in z. Let ifo = ^[1|2| ■ ■ ■ |n] be the weighted sum of CRSFs of Q (CRGs in 
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which each node is in a separate component). We have 



where au is the weighted sum of CRSFs with k cycles winding around the annulus. 

While we have not attempted to show that every connection probability can be computed 
via the .if -entries, we present here some case of small ni,n2. 

5.1 Annulus with (2,0) boundary nodes 

Suppose there are two nodes on one boundary and none on the other. Then ^i^2{z) counts 
connections to 1 from 2. There are only two topologically different configurations, which 



are illustrated in Figure 4 We have the following theorem. 





Figure 4: The two topologically distinct ways to connect the nodes on an annulus with (2, 0) 
boundary nodes. 

Theorem 5.1. ^\ogJ^i^2{z)\z=i is the probability that the LERW to 1 from 2 crosses the 
zipper. 

Proof. Let A^ , B^, and a^ (respectively) be the weighted sum of cycle-rooted groves which 
contain k cycles winding around the annulus, and in which nodes 1 and 2 are (respectively) 
connected by a path not crossing the zipper, connected by a path crossing the zipper, or are 
not connected. Then 

\2\ 



J^0 = Y,[Au{2 - z - l/zf + Bkzi2 - z - l/zf] = A^ + B,z + 0((z - 1) 

fc>0 

ir[l|2] = Y, «fc(2 - z - l/zf = a, + 0{{z - If) 



fc>0 



Bv [Theorem 4:21 ^1,2 = ^[?]/^[l|2], so 



^ \o^^ (z) - ^^^['] ^'^™ - ^° + o(z 1) n 

- log.i-i,2(^) - -^^ ^;p^ - -^^^^ + 0{z - 1). D 



5.2 Annulus with (1,1) boundary nodes 







Figure 5: Different numbers of windings of tlie patli to 1 from 2. These configurations 
contribute to Ci , Cq , c_i , and c_2 respectively. 

We let Cj denote the weighted sum of CRGs connecting 1 and 2 and such that the path 



from 1 to 2 crosses the zipper j times algebraically (see [Figure Sj ). Then ^[i] = ^j^iCj 



from which one can extract Cj for each j. ITheorem 4.21 shows that .ifi,2 = ■^[i]/^[l|2]: 
so =^1^2 and j!°[l|2] together determine the winding distribution (which we knew already 
from ITheorem 4."5l) . But j!°[l|2] is also a function of 2;, so =2i,2 does not by itself determine 
the winding distribution. But from .ifi^2 we can extract the expected number of algebraic 
crossings of the zipper via 

d 
E[# algebraic crossing of zipper] = — -log=Sfi 2 

oz 



5.3 Annlus with (2, 1) boundary nodes 

Suppose that nodes 1, 2 are on the inner boundary and 3 on the outer. Suppose the zipper 
starts between 1 and 2 (cclw) as in Figure [61 

This case can be reduced to the previous cases using [Theorem 4"!4l For example, consider 
the case that all three nodes are connected, and the path from 1 to 2 goes either left or 
right of the hole, as illustrated in [Figure 6[ 



Here by ITheorem 4.41 above. ^{^ is 1/^ times the sum of CRGs in which 1 and 2 are 
connected by a path, and 3 is considered an interior vertex (that is, either connected to this 
component or connected to a CRT disjoint from it); such configurations are weighted by the 
power of z determined by the path from 1 to 2. Thus ^i^3fo = A + Bz + C(2 — z — \/ z) 
and A and B are the desired quantities. 



?2,3 



jl,3 



The quantities ^{2 ^^"^ "^i'2 carry additional information about the winding of the 
path from 1 to 3 (respectively from 2 to 3) around the zipper. 
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Figure 6: The two main cases for connecting three boundary points on the annulus in the 
(2, 1) case (ignoring windings of the connection to vertex 3). 



5.4 Annulus with (3,0) boundary nodes 



This is a case which can be derived from the (2, 0) case using [Theorem 4.4l Suppose nodes 
1,2,3 are in counterclockwise order on the inner boundary, with a counterclockwise zipper 
between nodes 1 and 3. Consider the case when all nodes are connected; there are three 
possible configurations ^1,^12,^3 correspond to which complementary component of the 



2,3 



triple connection the outer boundary component lies (see Figure 7). The numerator of S^ ^ 3 



53,2 



is y4i + A2 + zA'i, the numerator of ££{2 is zA\ + ^2 + zA^^ 

zA\ + y42 + v43. These three quantities suffice to determine Ai,A2, A^ 



j3,l 



and the numerator of =2'2'i is 






Figure 7: The three topologically distinct subcases when the three nodes are connected (for 
the annulus with (3,0) nodes). 
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5.5 Annulus with (3, 1) boundary nodes 

This is a case we will need when we do the LERW computations 
below. Suppose there are 4 nodes in all, with nodes 1, 2, 3 on the 
inner boundary in counterclockwise order and 4 on the outer bound- 
ary. Suppose the zipper starts between 1 and 3 and is oriented 
counterclockwise, as in the figure. We wish to compute the ratios 
Z[l,2|3,4]/Z[l|2|3|4] and Z[l, 2|3,4]/Z[1, 2, 3,4] (as before, Z[(t] de- 
notes the weighted sum of groves of type a , for the trivial bundle) . 

Recall that ^[ila] denotes the weighted sum of cycle-rooted 
groves of type \\\ , times the parallel transport of the path to node 1 

from 2 and the path to node 3 from node 4, so that for a trivial bundle, ^[iH] = 2'[1, 2|3, 4] . 
Because of the path connecting nodes 3 and 4, in fact there will be no cycles in the CRG. 
Similarly, ^[l\i] and ^[Hf] denote the weighted sum of CRGs of type f|2 and 2I1 , times 
their respective parallel transports. 

Now by ITheorem 4.4t ^{2 has numerator counting connections f|| and ^[2 (with a 




?2,4 



minus sign). Similarly for =^^'3 and ^i'^- With J^ denoting the sum of CRSFs, we have 



^odet^^;^ 



=^ det .^2'i — ■^[2111 — "^[i 



ni\t] 
ni\t] 
lit] 



L3I1J 



^[2,4. 
1 
■3|4i 



nut] 



nut] 
nut] 



nut] 
ni 



2 1 4] 

sJ 



ni\t]-z'ni 



,312] 



(5) 



As a consequence 



nut] 



£e\ 



2.4 + ^2^M ^ ^3,4 



2,1 



When z — )■ 1 both the numerator and denominator of (jS]) converge to , so we evaluate the 
limiting ratio using rHopital's rule. We can expand J!fu,v = Lu,v + {z — l)-^u v + ^((^ ~ ^Y) 5 
where Lu,v is symmetric and L'^ ^ is antisymmetric. Then in the limit z — )■ 1 this gives 



Z[l,2|3,4] ,^^ nl\t] 



?2,4 I o (v?l,4 



3,4 



lim 



lim 



-^ 1,3 "T ^ -^ 3,2 "T "^ 2,1 



Z[l|2|3|4] z^i ^[l\2\3\A] z^i l-z^ 

-^1,2-^3,4 + -^1,2-^3,4 ~ -^1,4-^3,2 ~ -^1,4-^3,2 



-2 



+ 



+ 



(-^3,1-^2,4 + -^3,1-^2,4 — -^2,1-^3,4 — -^2,1-^3,4) + 2(1^3 iL2,4 — L2, 1-^3^4) 

^2 
-^2,3-^1,4 + -^2,3-^1,4 ~ -^1,3-^2,4 ~ -^1,3-^2,4 
^2 

"-^1,2-^3,4 ~ -^2,3-^1,4 ~ -^3,1-^2,4 + -^1,2-^3,4 ~ -^1,3-^2,4- 



(7) 
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It is useful to express this above formula in terms of the Green's function where the nth 

node is the boundary. We can express Lj^„ = — X]?=i -^ij ^^^1 iF[l, 2, 3,4]/iF[l|2|3|4] = 

det[^ij]ill'---'l^:l using ITheore^FTSl Let ^ij = dj + {z- l)G'ij + 0{{z - l)^) . Then using 

(|2]) to express =Sf in terms of ^ and taking the limit z — )■ 1, some algebraic manipulation 

yields 

Z[i,z|o,4j _ ^1 x-y/ n' \ n n fQ\ 

^M 2 3 41 ~ ~*^i,2 ~ '-^2,3 ~ '-^s.i + ^^1,2 - '^1,3- [o) 



(In Corollary 6.5 we will see a much easier way to convert an L-formula into a G-formula.) 
While the left-hand sides of these formulas ([7]) and (IHl) are symmetric in nodes 1 and 2 , 
the right-hand sides are not. This asymmetry is due to the location of the zipper, and moving 
the zipper would change the values of the L[j and the G'^j (albeit in a predictable way). 
If we keep the zipper between nodes 1 and 3, then we should expect a different formula 
for Z[l,3|2,4] than what we would get by permuting the indices 1,2,3 in the formula for 
Z[l, 2|3, 4] . Indeed, if we carry out the computations as above, we obtain 

— —L^ 2-^3,4 — -^2,3-^1,4 ~ -^3,1-^2,4 + -^1,2-^3,4 " -^1,3-^2,4 (9a) 



Z[l|2|3|4] 
Z[l,3|2,4] 



"-^1,2-^3,4 — -^2,3-^1,4 ~ -^3,1-^2,4 (9b) 



and 



Z[l|2|3|4] 

Z[2,3|l,4] _ , , , 

^r. |2|o|4l = "-^1,2-^3,4 - -^2,3-^1,4 " -^3,1-^2,4 + -^1,4-^2,3 - -^1,3-^2,4 (9c) 



-Gi2 — ^2,3 ~ G-^i + Gi^2 — Ci^s (10a) 

-g;,2 - G'2,3 - G'3,1 (lob) 

-G[,2 - G'2,3 - G'3,1 + G2,3 - Gi,3 (lOc) 



Z[l, 


2|3 


4] 


Z[l, 
Z[l, 


2,3 

3|2 


,4] 
4] 


Z[l, 
Z[2, 


2,3 

3|1 


,4] 
4] 



Z[l,2,3,4] 



5.6 Pair of pants with (2, 0, 0) boundary nodes 

The remaining two annulus cases are most easily viewed as special cases of the case when the 



surface S is a pair of pants with 2 nodes on one boundary and no other nodes, see [Figure 8 
Put an SL2(C) bundle with monodromies A and B around the two central holes Ca and 
Cb , and supported on zippers from the holes to the boundary between nodes 1 and 2 . 
The parallel transport of a path to node 1 from node 2 is of the form 

1. /, if the path has both holes on its right 
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A 

--{AB)0 A{AB)o 



B 

^{AB)OB{AB)" 



ABA-^ 

{AB)B{AB)- 



B-^AB ABAB-^A-^ 

=(AB)-iA(AB) =(AB)A(AB)-i 



Figure 8: Some of the possible topological types for the path between nodes 1 and 2 when 
the surface is a pair of pants with both nodes on one boundary. The lower zipper (in red) 
has parallel transport A, and the upper zipper (in blue) has parallel transport B. For each 
diagram, the parallel transport of the path to 1 from 2 is shown. 



2. AB , if the path has both holes on its left 

3. {AB)~'^ A{AB)'^ for some A; G Z, if the path has the lower hole on its left and the 
upper hole on its right, and k is the algebraic number of crossings that a dual path 
from the lower hole to the left boundary makes across the upper zipper 

4. {AB)^^B{AB)^ for some A; G Z, if the path has the lower hole on its right and the 
upper hole on its left, and k is the algebraic number of crossings that a dual path from 
the upper hole to the left boundary makes across the lower zipper 

We let c^^^-*, c^^^\ Cf, , and c^ (for A; G Z) denote the weighted sum of cycle- rooted 
groves of the above types. We further let c^ and c^ denote the number of cycle-rooted 
groves of type c*^^^-* and c^^^^ in which there are £ G N loops that surround both holes. 

We need to choose matrices A and B for which det A = 1 and det B = 1, and it is 
convenient to choose them so that Tt{A) = 2 and Tr(B) = 2 (so that loops which surround 
one hole but not the other have weight 0), and so that AB is diagonal. We can take 



A 



B 



2x 
x+1 



y 



2x 

X+l X 

{x-lfx 2 

j/(x+l)2 x+l 



2 

x+l 

y 



for variables x and y . Then we have 



AB 



X 





" 

l/x 
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and since AB is diagonal, it is straightforward to evaluate 



{AB)-''A{AB)'' 
{AB)-^B{ABf 



2x 
x+1 

y(x+lY 
2x 
x+1 

{x-lfx^+'^*' 
y{x+l)-2 



yx 



-2k' 



2 

x+1 

-yx~ 



-2fc-l' 



2 
x+l 



^[2] = /^cP(2 -X- l/xY + ABY,cf''\2 -X- l/xf 

+ ^cf''\AB)-'^A{ABf + ^cf''\AB)-^B{ABf. 



kez 



kei 



Only the last two sums contribute to the 1, 2 entry of J!°[f] : 



This is a Laurent series in x from which one can extract the coefficients c\, and c\. 
Once these are known, the coefficients c\ and c\ can be extracted from S'^n and 



5.7 Annulus with (2,2) boundary nodes 

On the annulus with 4 nodes, put nodes 1,2 on the outer boundary and 3,4 on the inner 
boundary. Suppose we wish to compute the probability of the connections 13|24 and 14|23. 
This computation can be used to compute the probability that an edge e is on the LERW 
from 1 to 2 (an equivalent calculation was done in [KenOOaj ). 

Insert an extra edge 634 from 3 to 4; this "splits" the inner boundary into two (See 
Figure [H]). This is then a special case of the construction of section EB In the notation of 
that section, it suffices to use x = —1 and x = 1 to distinguish crossings 13|24 and 14|23: 
[■^,2] 1.2 in the case a; = — 1 gives the sum and in the case x = 1 the difference of the two 
desired quantities. 

5.8 Annulus with (4,0) boundary nodes 

When 4 nodes are on the outer boundary and none on the inner (and nodes 1, 2, 3, 4 are in 
cclw order), the case we have not yet discussed is the 14|23 case: there are three subcases 
depending on whether the paths from 1 to 4 and 2 to 3 go left or right of the inner boundary. 
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Figure 9: Computing crossings 13|24 and 14|23. 

Again this is a special case of the (2, 0, 0) case if we put an extra edge between nodes 3 
and 4. 

In this case the only possible parallel transports to 1 from 2 are (using the connection 
from that section) 

ABJ,B,A,B-^AB, 

and only in the first two cases is there a possible extra loop surrounding both Ca and Cb ■ 
Thus 

^^^^ = I{c'i'''hcf'''\2-x~l/x))+AB{c^^''hcf'''\2-x-l/x))+Bc'i'''hAci^''hB-^ 

The 1, 2 entry in iFi_2 is 

[=^i,2ji,2 - --Co +yCo +-^q ■ 
From this we can extract the three cases of interest. 



6 Annular-one surface graphs 

Suppose n is even and that the first n — 1 nodes are on the inner boundary of the annulus 
arranged in counterclockwise order, and node n is by itself on the outer boundary, and that 
the zipper is between nodes n — 1 and 1 and directed in the counterclockwise direction (from 
1 to n — 1), as in section [5751 We call these annular-one surface graphs; they are the next case 
after circular planar graphs, and they play an important role in our loop-erased random walk 
calculations in section [HI Annular-one surface graphs of course include the (1, 1) and (3, 1) 
cases that we did in the last section, but for expository purposes we treated those special 
cases in separately. We are interested in computing, for any partition a, the weighted sum 
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of groves of type a, which we denote Z[a]. We show how to compute Z[a]/Z[l\2\ ■ ■ ■ |^] in 
terms of the response matrix ^, and Z[a]/Z[l,2, . . . ,n] in terms of the Green's function 
^ . We treat first the case when a is a complete pairing, and then show how to reduce the 
computation for general partitions to the computation for pairings. 

6.1 Complete pairings 

There are n — 1 ways to connect the two boundaries (ignoring windings). When the annulus 
is cut along this connection, the domain becomes planar, so there are Cn/2-1 ways to pair 
up the remaining nodes. We have 

{n - l)C„/2-i -{n- 1)^^/2 -l)!(n/2)! ~ 2 (n/2)!(n/2)! ~ 2 [n/2 

So the number of annular pairings equals the number of equations arising from the deter- 
minant formula. In fact, there is a natural bijection between the ^-determinants and the 
annular pairings which is based on the cycle lemma of Dvoretzky and Motzkin |DM47] . 
which we use in Appendix A.l to show that these equations are linearly independent for any 



even n. 

The approach to computing the normalized probability of pairings is similar to the (3,1) 
case above. In any directed pairing, the connection between node n and the other boundary 
determines whether or not and in what direction that any other directed pair crosses the 
zipper. Reversing the direction of any directed pair (other than the pair containing n) 
that crosses the zipper introduces a factor of z^ . Since only even powers of z appear, it is 
convenient to change variables to 

For example, ^[l\l\l] = C^[il?l2- 

For more compact notation let ^o- := j^rmf |„i ■ When expanding an =Sf -determinant 
into a signed sum of ^o- 's, where a is a directed pairing, ^„ can be put into a canonical form 
^power^^^^ where a' is a directed pairing in which the pairs are directed counterclockwise 
around the annulus. We can then solve for any given ^r = S't{z) in terms of the =Sf- 
determinants and C, and take the limit C, ^ 1. 

The system of linear equations can be represented by a matrix A„ . When recording the 
linear equation corresponding to det =Sf^ , we can re-order R and 5* in any manner, and this 
would just scale row det^SfJ^ of A„ by ±1, which has no effect on our ability to solve for 
the ^[cr]s. But the signs in A~^ are surprisingly nice when we order R and S* in a manner 
that corresponds to det=Sf^'s associated pairing in the aforementioned bijection. 

The matrix for A2 is just the 1x1 matrix whose entry is 1 (since ^[^]/^[l|2] = ^12 = 
det^/): 
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det ££\ 1 
The matrix A4 encodes the system of equations dS]) we saw for the (3, 1) case: 



;^ 



det^i'3 



det^. 



1,4 
3,2 



det=^^'^ 





COtM 





-1 


1 





-c 


1 



The next matrix Ar is 



coin 


^10 


to-^ 


to^ 


;OC0 


^CO 


cocq 


cooq 


to^ 


to^ 


^00 


CO(M 


CO(M 


(M^ 


Cqrt 


^\a 


^10 


10^ 


lO-!* 


^co 


cq^ 


^rH 


^LO 


CO 10 


10^ 


O)-* 


'^CO 


^CO 


cocq 


incq 



;^ ;^ 



det^^'^'^ 

ueL,=z- 1,3,5 
det =£ \^2,'o 
det =^5'2|4 

Cie u oZ- c; 1 A 

oe u rJZy 4 1 Q 

det^^'^'^ 
ueL=i 453 

det =z:^^^2 

det =-2^3'4|2 
Oe u c^ 2 4 1 

det =^2!3',i 



1 


-1 








1 





1 


-1 


1 


-1 














c 





1 


— 


-1 























— ( 















-1 
1 

-C 



C 





-1 



1 

1 

-1 





















-1 



1 


c 
-c 





c 













1 

-1 



-1 



1 




-1 


1 





-1 
1 



Each det=^J^ is a signed-sum of {n/2)\ of the iF[o"]'s, but not all of these a's can be 
embedded in the annulus, so the rows generally have fewer than {n/2)\ nonzero entries. 
For each pairing a that embeds in the annulus, the column ^[cr] contains 2"/^~^ nonzero 
entries: for each pair {2,7} in a, except the pair containing n, either i & R and j G S* or 
else j E R and i E S . 
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The inverses A„ ^ of these matrices are 



^1?] 






^ CO ^ (M '^ ^ 

cs" ^ r-T CO CO c4^ 



3I2J 



lit] 



w[ 



cu 






1 c 
1 1 

c c 



cu 



X 



;i-c)^ 





to 10 
^[^co^ 

cu 


CO 10 
CO (m" 

OJ 

-a 


CO cnT 

cu 


CO 10" 

cu 


CO CO 

cu 


CO CO 

cu 


CO O) 

^"co" 

cu 


cu 


CO ^ 

co'^of 

cu 

1I3 


CO rt 

^"co" 
in" of 

cu 


■^■[2 4161 
-=^ [1 3I5J 


"C+1 


c+1 


C^+C 


2C 


c+1 


c^+c 


2C 


2C 


c+1 


2 


•^■[4 3161 

■=<f Li 2I5J 


c 


1 


c 


c 


c 


c^ 


c 


c 


1 


1 


■^■[1 3161 
=<f [5 2I4J 


c+1 


2 


c+1 


c+1 


c+1 


2C 


c+1 


c+1 


2 


2 


t^lS 2 61 

=<f [5 1 4J 


1 


1 


c 


1 


1 


c 


c 


c 


1 


1 


■^■[5 2161 

=^U 1I3J 


2C 


2C 


c^+c 


2C 


c+1 


c^+c 


c^+c 


2C 


c+1 


c+1 


tfk=''r2 1 1 61 
=<f [4 5I3J 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


■^■[4 1161 
=^ [3 5I2J 


C+1 


C+1 


2C 


2C 


c+1 


2C 


c+1 


c+1 


c+1 


2 


'(W'\i 5 161 
^[3 4I2J 


c 


c 


c 


c 


1 


c 


c 


1 


1 


1 


■^■[3 5161 
-^ [2 4I1J 


c^+c 


2C 


c^+c 


c^+c 


2C 


2e 


c^+c 


2C 


c+1 


c+1 


'tfv'rs 4 161 
■^[2 3I1J 


c 


c 


c^ 


c^ 


c 


e 


c 


c 


c 


1 



X 



(1-C)^ 
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J^[l|2|3|4|5|6] (1-C)' 



7l,3,6 



For example, the fifth row tells us 

2C det if ?; J^ + 2C det if fj' + (C' + C) det ^^5,2,4 

+2Cdetifg^ + (C + l)deti^g^ 

+ (C^ + C) det if J^'l + (C^ + C) det if^'j;^^ 

+2C det i^^ J^ + (C + 1) det if |g + (C + 1) det i^^;? 



) (11) 



and the sixth row tells us 



n 



2|l|6l 
4I5I3J 



det i^it^ + det ^t'2'5 + det i^5'2'4 



^[1|2|3|4|5|6] (1-C)^ 



+ det ^t:TA + det ^TA + det ^YA + det ^: 



4,5,3 



3,5,2 



:i2i 



5,4,6 
2,3,1 



+ det ^^ JJ + det =^m + det ^ 
Theorem 6.1. For positive even integers n, the matrix An is nonsingular. 

We postpone the proof of this theorem until the appendix. 

There are some clear patterns in A~^ which warrant further investigation: 

Conjecture 6.1. For positive even integers n, the entries of the matrix (1— C)"/^~^x A^^ are 
all polynomials in C, , with nonnegative integer coefficients, and with degree at most njl — 1 . 
These polynomials, when evaluated at ( = 1 , only depend on the row of the matrix, and 
always divide (n/2 — 1)!. These numbers count "Dyck tilings" (certain objects defined in 
IKWllbf ) whose lower boundary is the Dyck path associated with the pairing that indexes 
the row of A^ ■ (We have verified this for all n < 14.) 

Conjecture 6.2. For positive even integers n. 



det A. 



:i-c) 



on — 2_ _L 



2 V'i/27 



(We have verified this for all n < 16 .j 

Proposition 6.2. For positive even integers n, the row in (1 — (^)"/^^^A^^ corresponding to 



■^■rn/2-l|n/2-2| I 1 I „ 1 

•^ Ln/2+1 In/2+2 I ' ' ' |n-l|n/2j 



consists of all 1 's. 



Proof. We compute the all ones horizontal vector multiplied by A„. The resulting vector 
is indexed by pairings a . If the matrix entry {An)(A,B),u is nonzero, then a pairs A to B . 
If a contains a pair which does not cross the zipper, then let (i, j) be the lexicographically 
smallest such pair. We can "swap" the pair (z, j) and find that 

(An)((A\{i,j})U{i},(iJ\{i,i})U{j}),a = -(A„)((A\{i,i})U{j},(iJ\{i,i})U«),a, 
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so the a entry of the product vector is zero. If every pair of a crosses the zipper, then a 
must be do = ^^2+1 In/2+2 1 ■ • • \n-i\n/2- We havc 

(A„)((A\{ij})u{i},(i?\{i,j})u{j}),ao = -C X (A„)((A\{ij})u{i},(B\{i,j})u{i}),<7o, 
for i < j , and hence the product vector is (1 — C)"^^^^ at (Jq. D 

We can expand out the ^-determinants into sums of products of the =Sfj.j , each of which 
depends on (, where Jfj^i{C) = J^ij{l/C). Since the denominator is (1 — C)"^^^^; to evaluate 
the hmit C ~^ 1? we can differentiate the numerator and denominator n/2 — 1 times with 
respect to ( and then set C to 1. The denominator of course becomes (— l)"/^~^(r2/2 — !)!• 
The numerator will consist of monomials of degree n/2 in the quantities 



^- I -if- ■ 






C=i 



nn/2-1 
Q(-n/2-l'^''-'3 



C=l 



In each case all the terms involving higher order derivatives of =Sfjj- cancel upon setting C, 
to 1 . This is quite convenient, since there are fewer quantities that we need to evaluate later. 
We already have 



Lid 



^ -I 



^- -I 



let us define 



L' 



«j 



dz •■' 



z=l 



d_ 



'«J 



c=i 



Theorem 6.3. For all positive even n and pairings o o/ {1, 
nomial of degree n/2 in the quantities 



n], Z„/Zi\2\-\n IS a poly- 



{Lij : I < i < j < n} and {L[j : 1 < i < j < n — 1}. 
We prove this theorem in the appendix. 



Conjecture 6.3. The coefficients in the polynomials in \Theorem 6.3\ are all integers. (We 
have verified this for all a for all n < 10.) 
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In the above two examples we obtain 

^[4.5|1,2|3,6] ///+///+// L 

Z[l|2|3|4|5|oJ 

— -^1,2-^3,5-^4,6 ~ -^3,6-^4,5-^1,2 + -^3,5-^4,6-^1,2 + -^2,6-^4,5-^1,3 " -^2,5-^4,6-^1,3 

~ -^2,6-^3,5-^1,4 + -^2,5-^3,6-^1,4 ~ -^1,6-^4,5-^2,3 + 2Li,5L4,eL2 3 — Li/^L^^qL2'^ 

+ 2-^5,6-^^4^23 ~ 2^4,6-^15-^2,3 + -^1,6-^3,5-^2,4 ~ 2^1,5^3,6-^2 4 + -^1,3-^5,6-^2,4 

— 2L5,6-Z^i 3-^2,4 + 2-^3,6-^1,5-^2,4 + -^1,4-^3,6-^2,5 ~ -^1,3-^4,6-^2,5 + 2-^^4,6-^1,3-^2,5 
~ 2L3,6-Z^1 4-^2 5 — -^1,6-^2,5-^3,4 + 2-^1,5-^2,6-^3,4 ~ -^1,2-^5,6-^3,4 + 2-^^5,6-^1,2-^3,4 

— 2L2,6-Z^i,5-Z^3,4 + 2-^1,6-^2,5-^3,4 ~ -^1,4-^2,6-^3,5 + -^1,2-^4,6-^3,5 ~ 2-^4,6-^1 2-^3,5 
+ 2-^2,6-^1,4-^3,5 ^ 2-Li,6-L2,4-^3,5 + -^1,3-^2,6-^4,5 ~ -^1,2-^3,6-^4,5 + 2-^3,6-^1,2-^4,5 

— 2-^2,6-^1,3-^4,5 + 2-Li,6-/^2,3-^4,5 

-^5,6-^1,4-^2,3 ~ -^4,6-^1,5-^2,3 ~ -^5,6-^1,3-^2,4 + -^3,6-^1,5-^2,4 + -^4,6-^1,3-^2,5 
^-^3,6-^1,4-^2,5 + -^5,6-^1,2-^3,4 ~ -^2,6-^1,5-^3,4 + -^1,6-^2,5-^3,4 ~ -^4,6-^1,2-^3,5 
+-^2,6-^1,4-^3,5 ■" -^1,6-^2,4-^3,5 + -^3,6-^1,2-^4,5 ~ -^2,6-^1,3-^4,5 + -^1,6-^2,3-^4,5 



Z[4,2|5,l|3,6] 
Z[l|2|3|4|5|6] 



This last formula can be written as a Pfaffian: 

-^1,6 
-^2,6 
-^3,6 
-^4,6 
-^5,6 
_ 

Notice that the derivatives appear everywhere except the last row and column. 










T' 

-'^1,2 


-^1,3 


^'1,4 


^1,5 






-^^1,2 





Tl 

^2,Z 


^'2,4 


-'^2,5 


Z[4,2|5, 13,6] 


- Pf 




Tl 





Tl 

-^3,4 


Tl 

-^3,5 


Z123|45 6 




-M,4 


-^2.4 


-'^3,4 





-^4,5 






-A,5 


-^2,5 


_TI 

-^3,5 


-^^5 









~-^l,6 


~-^2,6 


—-^3,6 


~-^4,6 


-^5, 



6.2 Partial pairings 

When we carry out our loop-erased random walk computations, we also need to handle 
partial pairings, where the grove partition has singleton components as well as pairs, so 
long as node n is in a pair. The determinant formulas from [Theorem 4.41 are, as formal 
polynomials in the matrix entries of the response matrix ^ , exactly the same as they would 
be if the singleton parts of a were removed, and these singleton nodes were made internal 
vertices. As a result, the "-L-polynomials" look the same, even though the underlying L- 
variables depend on the set of nodes. For example, to find Z[l, 2|3, 5|4] we can take fl9aj) 
and simply replace the 4-indices with 5's: 

Z[l,2|3,5|4] 



Z[l|2|3|4|5] 



"-^1,2-^3,5 — -^2,3-^1,5 — -^3,1-^2,5 + -^1,2-^3,5 — -^1,3-^2,5 
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The formula in terms of G and G' does however look different. 

Theorem 6.4. Let a he a partial pairing, in which the nodes T are in singleton parts, with 
n ^ T . In the above formulas expressing ^[o"]/iF[l|2| ■ ■ • |^] in terms of ^ -determinants, 
we can replace each det ^^ with det ^J'y J , where '^ij = '^ij for j < n and ^i^n = 1 ; and 
the result will be a formula for Z[(t]/Z[1,2, . . . ,n]. 

Proof. Since the higher derivatives of Lij and the first derivative of Lj„ always cancel out, 
we may compute Z [o"] using any convenient choice of L^ „ , and for present purposes it is 
convenient to make the choice for which ^ . ^ij = for each i. Then writing the sequence 
S as S = S*,n, where n ^ S* , we may express the determinant det=SfJ^ as 

det^J' = -^det^J'\ 

Let T be the sequence {1, . . . , n — 1} \ (i? fl S**) in sorted order. Suppose for now that R 
and S are also in sorted order, and that R,T and S*,T are also in sorted order. We can 
use Jacobi's formula on each summand to obtain 

^^^■^R _ _Y^(_.^J2R+Y:S*+^+\{s(^S'■■.s>^}\ , .( cfi-l-^{l,...,n-l}\R 

Jpt 9^1-..,n-l - Z^^ ^) ^^^^-^ ){l,...,n-l}\{S',ir 

Now we use {1, . . . ,n — 1} \ R = S*,T and {1, . . . ,n — 1} \ {S*,i} = R,T \ {i} , and the 
fact that .^i '.[['nil and ^^^ [',','n~i are negative inverses to write 

det^l 






l,...,n-l V ^J /_^ V "-) ^^'"^R,T\{i}- 

ieRUT 



When we expand detSf^^""' along column n, we obtain 



\RUT\ 

' T\j+\R\ A^i-CH 

R,T with jth item removed 



det^l*/'^ = J2 (-l)^+l''l det^^^*'^ 



If the j'th item oi R,T is i, then i — j = \{s E S* : s < i}\, so 

\RUT\ 

det^l*/'^ = (-l)l«l+l^*l J2 i-iy^^^'^'*--'^'^Uet^[^l^^^ 



and so 



det^J__ ^ f^^j:R+j:S' + \R\ + \T\ + \R\ + \S*\^^^c^S*,n,T 



det^;:::;::^~' ' ^-^ ' 
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and since RU S* = {1, 2, . . . , 2|i?| — 1}, which adds up to \R\ modulo 2, we have 



"^""^^^ det^^'^. (13) 



det -^1 



l,...,n-l 
,n-l 



Now observe that if we relax the assumption that R, S , R,T , and 5**, T are in sorted order, 

the left- and right-hand sides of the above equation change signs the same number of times. 

Finally, recall that ^L^r":" = det —.^i^',',',]n-i ■ D 

Corollary 6.5. For a complete pairing a , the Green's function formula for Z[o]/Z\i, 2, . . . , n] 
can he obtained from the response-matrix formula for Z[o]/Z\l.\2\ ■ ■ ■\n] simply by replacing 
each Lij with Gij and each L'^ • with G'^ • , and then setting Gi^n = 1 ■ 

Corollary 6.6. For a partial pairing a in which node n is paired, the Green's function 
formulas for Z[a]/Z[l,2, . . . ,n] are invariant under the addition of global constant to the 
Green's function. 

Proof. The column indexed by n in def^^'j, is all-ones. D 

6.3 Windings 

We can also extract information about the windings of the paths within a grove pairing in a 
manner similar to that described in the (1, 1) case. For a given directed pairing a, we have 

k 

where Z[o", (/c)] is the weighted sum of groves of type a in which the algebraic number of 
zipper crossings (involving all pairs in a) is k. Then 

d 2f\a] 
E[algebraic number of zipper crossings for groves of type a] = lim — — log -— . 

z^i dz J^[l|2| ■ ■ ■\n\ 

6.4 General partitions 

For more general partitions o", we can express Z[a\ as a linear combination of Z[t\ 's where 
the r's are partial pairings which may have unlisted nodes. We explain first how to do this 
with circular planar graphs (where the surface is a disk), and then for annular-one graphs. 
It is also possible to do this reduction for the annulus with 2 nodes on each boundary, but 
we do not know how to do it for the annulus with 2 nodes on one boundary and 3 on the 
other. 

Consider the disk with nodes labeled l,...,n in counterclockwise order on the outer 
boundary. For a partition a, let i be the smallest node label that is in a part of a of size 
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more than 2, and let s be the size of this part. We will measure the "badness" of the 
partition a by the quantity n{n — i) + s. Let j be the next-smallest item in i's part of a. 
Let a* denote the partition obtained from a by "de-listing" j, i.e., by regarding j as an 
internal vertex which can occur in any of a* 's parts. If a has k parts, then we can write 
Z[a*] = ^£=1 Z[a* with j added to ith part] . One of these terms on the right will be Z[a] , 
so 

Z[a] = Z[a*] -^Z[a* with j added to ith part]. 
e 

where the sum runs over all parts of a* except the one containing i. Because the graph is 
circular planar, unless j is added to a part of a* that is "covered" by the part containing i , 
there will be no groves of that partition type. Thus each of the nonzero terms on the right 
has a smaller badness than a, so we can iterate this process to eventually express Z[a] as a 
linear combination of Z[t] 's where r is a partial pairing. For example, 

Z[l,5,8|2,3,4|6,7] = Z[l, 8|2, 3,4|6, 7] - Z[l, 8|2, 3,4, 5|6, 7] - Z[l, 8|2, 3,4|5, 6, 7] 
= Z[1,8|2,4|6,7]-Z[1,8|2,5|6,7]-Z[1,8|2,4|5,7]. 

For an annular-one graph, we can do essentially the same procedure as for a circular 
planar graph, except that we start with the part containing node n, reducing its size if the 
part has more than two nodes. Once this part has size two, say that it is {h,n}, then we 
list the remaining nodes in the order h + l,h + 2, ... ,n — 1,1,2, ... ,h — l, and do reductions 
described above. If any node gets adjoined to the part {h, n} , then the term for that partition 
will be zero. If n started out in a singleton part, we can start out as in the circular planar 
case (with the order 1, . . . ,n— 1) until a node gets adjoined to n's part, and then cyclically 
re-order the nodes as in the doubleton case. 

7 The Green's function and its monodromy-derivative 

To carry out our loop-erased random walk computations for various lattices, we will use our 
formulas for the connection probabilities in annular-one graphs developed in section El and 
for this we need the Green's function G together with its derivative G' with respect to a 
zipper monodromy. We will need G and G' for both the full lattice, and the lattice after 
some of its edges have been cut. 

7.1 Green's function and potential kernel 

The Green's function Gu,v is infinite for recurrent lattices such as 7? , but there is a quantity 
known as the potential kernel Au^v which behaves like a Green's function, except that Au,u = 
0, and Gu,v and Au^^ have the opposite sign convention (see |Spi76] ). Suppose that a graph Q 
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is the intersection of Tl? or another lattice L with a region surrounding the origin, with 
"wired boundary conditions", i.e., all the lattice vertices in L that are not in the region are 
merged into a single vertex in Q that plays the role of boundary. If R denotes the electrical 
resistance within Q from the origin to the boundary, then 

Gj, = i?-<, + o(l), 

where the error term tends to for fixed u and f as i? — > oo. For translation invariant 
lattices, J^ ^ depends only on u — v, and is written as a^_^ . 

Since all of our formulas for crossings of the annulus are invariant when a global constant 
is added to the Green's function (involving terms such as Gi^2 — Gi^s), it is straightfoward to 
take the limit hmg_>.L of these formulas by replacing each G^^ in the formula with —A^^, 
which we shall also denote by Gj;^,. 

For convenience let us work with a modified finite graph Q approximating the lattice L, 
obtained by adjoining an edge with conductance —1/R to node n oi Q , taking node n of 
Q to be the other endpoint of this edge. This has the effect of making the resistance in Q 
from the origin to node n exactly zero. 

For any partition a for which n is not in a singleton part, we have Z^ = —Z/R and 
Z^[a] = —Z[a]/R, so in particular we can compute Z[a]/Z = Z^[a]/Z^ by working with 
the Green's function ^^ of this modified graph. We define Gu,v and G'^ ^ by the expansion 

<. = G.,. + (^-l)G:,, + 0((z-in. 

Then in the limit ^ — )■ L, we have Gu,v -^ ~^uv- 

It is well-known how to compute the potential kernel on periodic lattices by taking the 
Fourier coefficients of the characteristic polynomial of the lattice |Spi76| . The potential 
kernel can also be computed for any "isoradial" graph by doing local computations |Ken02j . 
The square, triangular, and hexagonal lattices are both periodic and isoradial, so for these 
lattices either method can be employed. 

We shall make use of the following smoothness result: 

Lemma 7.1 ( |Sto50t lKen02j ). For points z = (21,22) far from (0,0), the potential kernel 
on T? behaves like 

ZTT /TT 

where 7 is Euler's constant. 

Similar formulas hold for other lattices |Ken02] . 
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7.2 Derivative of the Green's function 

7.2.1 Infinite sum formula 

Let S be the adjacency matrix of t lie zipper, i.e., 

J 1 there is a zipper edge directed from k to i 
I otherwise. 

Then A{z) = Aq + (1 - z-^)S +{l-z)S*, so 

A{z)-' = (Ao(l + (1 - z-')A,'S + (1 - z)A,'S*))-' 

= A,' - (1 - z-')A,'SA,' - (1 - z)A,'S*A,' + 0((^ - 1)^) 

'^u,v = Gu,v ^ (^ ^ 1) 2_^ Ck,e{Gu,kGe,v — Gu/Gk,v) + 0{{z — 1) ) 

zipper edges {k,£) 

The sum is over zipper edges {k,£) in which the zipper direction is from k to i, and Ck/ is 
the conductance of edge {k, i) . The linear term in z — 1 gives us the desired derivative: 

Gu,v — ^z^u,v\z=i = ~ 2^ Ck/{Gu,kGe^v — Gu/Gk,v)- (14) 

zipper edges {k,l) 



For the modified graph ^ , we of course have 



^U,V ~ ^Z^UM 



2=1 

zipper edges {k,l.) 



7 ^ Ck/{Gu,kGi^v — Gu/Gk,v)- (15) 



For a vertical zipper in 7? (or the triangular lattice or hexagonal lattice) started in the 
face whose lower-left corner is the origin, directed down towards infinity, we define 

G'u,v = ~ 2-^ '^k,iiGu,kGi^v~ Gu/Gi^,^^). (16) 

zipper edges (k,£) 

For fixed u and v, for zipper edges {k, i) at a distance r from the origin, it is straightforward 
to use the smoothness result in lLemma 7. II to show that edge {k,i) contributes 0(r~^logr) 
to the sum[T6| so this sum is absolutely convergent. 

We would like to know that G'^^, converges to G'J;^, as defined in flT6l) for a sequence 
of ^'s converging to L. For our purposes in section [8] when we analyze loop-erased random 
on the lattice, we do not need this convergence of G' for every sequence of ^'s converging 
to L, it will sufiice to have convergence for some sequence of ^'s tending to L. Perhaps 
the easiest way to show this is to exploit the refiection symmetry that each of the square, 
triangular, and hexagonal lattices possess. 
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Lemma 7.2. If h is the square, triangular, or hexagonal lattice, and L G N, let Ql = 

[—L^, L^] X [—L, L'^]nL be the off-center box surrounding the origin and zipper (as in Figure 3 



except with the lower boundary much closer to the origin than the other boundaries) , where 
the lower boundary of the box is aligned with an axis of reflection symmetry of the lattice L. 
Let u, V be fixed points in L. Then 



lim G J^^ — G ^ . 

Proof. We can approximate G^^ and G^^ (for w within distance L of the origin) using 
the Green's function of the lattice intersected with the upper-half plane. More precisely, we 
approximate G^^ by 

G~ •= —A^ 4- A^ 

p,q ' P^g p* ,q^ 

where p* = p — (0, 2L) is the reflection of p through the lower boundary of the box, and A^ 
is the potential kernel of the lattice. By construction G^ ^ is zero for q along the lower side 
of the box, and by the smoothness result from lLemma 7.ll Gp = 0(1/ L"^) along the other 
three sides of the box. Both G^ ^ and G^^ are harmonic in both p and q within the box 
(except on the boundary), and G^^ is zero along all four sides of the box. By the maximal 
principle for harmonic functions, for p and q within Ql we have 

\Gp^g - Gp^^l = 0{1/L ), 

i.e., 

Gp^q = ~\,q + \*,q ~ ^0*,0 + G{l/L ). 

Next we compare the contribution of a zipper edge {k,i) to G'f^% and G'J;^. We have 
~\^u,k^e,v ~- ^u,e^k,v) 



— f—A'^ 4- A^ — A^ 


o)( "^v/ + ^v 


- A^ 


o) 






_|_ ('—/iL _|_ ^L _ aL 


){-Al, + Al 


- A^ 

',k ^0* 


o) + 


0{L~Hog 


L) 
















.,J 










1 /iL /iL _|_ aL Ah _ 


_ /iL /iL _ 
^u/^v*,k 


Ah Ah 

^u*,e^v,k 








_i_ /iL r 4^ _|_ /iL 


Ah 1 Ah 


^ A^ - 


At 


-^ A^ - 


Al 


+ OiL-HogL) 












- ['~'u,k'~'e,v - ^ui^k.v) 


( nih (^h 

y^u,k*^f*,v 


-Gl,X 


ih \ 

k*,v) 







I /iL (Ah _ Ah \ _ ( Ah _ Ah \Ah 



Ah (Ah _ /iL \ . (Ah _ /|L \Ah 



+ 0{L-HogL) 
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Recall that u and v are fixed, so they are within distance 0(1) of the origin. The second 
term is 0{L~^ log L) . If the zipper edge {k,i) is at distance r from the origin, then the 
next four terms largely cancel one another and add up to 0{l/{rL)) . Upon summing over 
all zipper edges, we find G'% = G%^ + 0{L-^ log L) . D 

7.2.2 Zipper deformations 

The next task we have is to evaluate in closed form the infinite sum in fll6p . This we can do 
for many lattices L, including the square lattice, triangular lattice, and hexagonal lattice, 
although it is not clear how to do this for arbitrary lattices. 

We shall need to deform the path that the zipper takes. In general deforming the zipper 
while keeping its endpoints fixed has no effect on G'^^, unless the zipper is deformed across 
either u or v. If the zipper is moved across u in the direction of the arrow on the zipper, 
then G'^ ^ decreases by Gu,v , and similarly, moving the zipper across v (in the direction of 
the arrow) increases G'^^ by Gu,v ■ We can also move the endpoint of the zipper by adding 
a new zipper edge {k, i) (or removing an old one) near the endpoint of the zipper, which of 
course just adds (or removes) one term to the summations 0151) and 0161) . 

7.2.3 Closed- form evaluation of G' on Z^ 

Next we evaluate G' for Z^. The first step is to rotate the entire lattice 180° about the 
terminal square of the zipper. The rotation of course preserves the lattice Z^ , and maps u 
and V to (1, 1) —u and (1, 1) — v respectively, but now the zipper goes up to infinity rather 
than down to infinity. Let G'^ denote G' with the repositioned zipper. We have 

'^U,V — '-^(1,1)-U,(1,1)-!,- 

The next step is to deform the zipper so that it once again goes downwards. We can 
deform the initial segment of the zipper so that it goes downwards, then circles around back 
up along a large-radius circle, and then continues back up as before. By ILemma 7.21 the 
summands along the zipper starting with the large-radius circle and the subsequent path to 
infinity are negligible. So we have 

^'u,v = G'li,iyu,{i,i)-v = G[i^i)-u,{i,i)-v + another term, 

where "another term" refers to the term for the zipper being deformed across u or v . 

Next we move the location of the start of the zipper, translating it by v + u — (1,1), by 
adding a finite number of new zipper edges. Then we deform the zipper again, making it go 
straight down; we have to add another term if the zipper gets deformed across either (1, 1) —u 
or (1,1) — V. Because the lattice Z^ is invariant under such translations, translating the 
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starting face of the zipper is equivalent to translating the vertices in the opposite direction. 
Thus we have 

^'u,v = ^'^,u + finite number of terms. 

Finally we use the antisymmetry of G'^^: 

-I finite number of terms 

This procedure is perhaps better explained by way of an example. We can write 

'-^(0,0), (2,1) — "-^(1,1), (-1,0) 

= '-^(1,1), (-1,0) + G'(i,i),(-i,o) 

= '-^(2,1), (0,0) + (^(1,1), (0,0)^(0,1), (-1,0) - (^(1,1), (0,1)^(0,0), (-1,0)) + G'(i,i),(_i,o) 

_ (<J(i,i),(o,o)G'(o,i),(-i,o) - <^(i,i),(o,i)G'(o,o),(-i,o)) + G'(i,i),(-i,' 



fii _ ii') _ i _|_ 2 

VyrTT 4 4/ 4~'~7r 

2 
_ 1 1 _ 5 

~ 2^^^^ 32 



In like manner we can compute G"„ ^ for any pair of vertices u and f in Z^ . The 
will always be in Q + ^Q + ;^Q. 



answer 



7.2.4 Closed-form evaluation of G' on the triangular lattice 

We can compute G' on the triangular lattice in essentially the same manner as for 1? . The 
key properties of the lattice that we used is that it is invariant under 180° rotations, and 
that for any pair of vertices there is a lattice-invariant translation that maps the first vertex 
to the second vertex. 

7.2.5 Closed-form evaluation of G' on the hexagonal lattice 

The hexagonal lattice is invariant under 180° rotations and is vertex-transitive. However, 
there are not lattice-invariant translations between any pair of vertices: we can partition the 
vertices into two color classes, black and white, such that any lattice-preserving translation 
will map the black vertices to the black vertices and the white vertices to the white vertices. 
Suppose M is a black vertex and f is a white vertex. After a 180° rotation about a 
hexagon, u is mapped to a white vertex m' and v is mapped to a black vertex v' . We can 
then translate m' to v and v' to u while preserving the lattice. This allows us to compute 
G'^ y when u is black and v is white (or vice versa) . 
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Since G'^^ is harmonic in both u and v (except along the zipper), when u and v have 
the same color, G"„ ^ can be expressed as 

r" = —(c' -^ c" A- c" ") 

^-'u,v q V U,U)1 "•" U,W2 u,wz) 

(plus another term if one of the edges {v,Wi) crosses the zipper), where the w^'s are the 
neighbors of v , and the right-hand side we can compute by the above method. 

7.3 Cutting edges 

Suppose that in the vector bundle setting, we know the Green's function ^ = ^^ for a graph 
Q , and we wish to know the Green's function for the graph G\{s., t} obtained by deleting an 
edge {s,t} of Q. Recall that Cs^t denotes the conductance of edge (s,t), and let us denote 
by r the parallel transport to s from t, so that Af^ = —Cg^t'T and Af^ = —Cs^t'T* ■ Then it 
is readily checked that 



a 



s,t 



where 



a 



s,t 



%,s + %,t~%,tr*-r%s-l/cs,t 



(which is a scalar). Indeed, if we let f{u,v) denote the purported Green's function on the 
right-hand side of (ITTl) . then f{u,v) = when either m or t> is the boundary, and we have 

J\u,v)/\^ = 6u,yj ■ . 

li w y^ s and w ^ t then Av}w' = A^^, so 

J2 /(«, ^)A^}j^'*^ = 5.,^ (iiw^s and wj^t). 
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Suppose now w = s. Then 



(^u.s - 




\^^ {%,s-%,tT*Ws,t- 


-r%,ty 


ct,sr* 


Ols,t 






— 




-^^m)1 


Ct,s 


(^n.s- 


Ct,s 
Ois,t 




- %,s + r%,s) 



by the choice of as,t ■ The case w = t is similar. 

Let us return to the hne bundle setting, with a zipper monodromy of z, that we are 
interested in near z = 1. If (s,t) is a zipper edge then r = 2; or r = 1/z (depending on 
the zipper direction), and otherwise r = 1. Let r' = dzT\z=i. Recall that G = ^\z=i and 
G' = dz^\z=i- From (^ it is evident that 



/^g\{s,t} _ ^ _ {Gu,s — Gu,t){Gs,v — Gt^i 



(19) 



''S,t 



where 



We have 



0-3,1 — Gs,s + Gt,t ~ 2Crs,t ^ l/'^s,t- 



dzas,t = dz'^s,s + dz%,t - idz%,t)r* - ^A^* - {dzr)%s - rdz^t,. 



(20) 



dz(ys,t\z=i = 0- 

Using this, we can differentiate (fT7|l with respect to the zipper monodromy z and set z = 1 
to obtain 



(nG\{s,t}y ^ r;' 



{^u,s ~ ^u,t + Gu,tT ){Gs,v — Gt,v) + {Gu,s — Gu,t){Gg,^ — G^^^ — r Gt,v) 



a-si 



One final edge-cutting formula is 

^G\{s,t} 



1 ~ Cs,t\Gs^s + Gt,t — '^Gs,t)- 



(21) 



(22) 
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This holds because directed edge (s,t) occurs in a uniform spanning tree of Q with proba- 
bihty Cs^t{Gs^s — G^s,t), and hkewise directed edge (t, s) may occur in the spanning tree. 

These formulas (fT9l) . (J20l) . (J2T1) and (J22l) of course apply to the modified graph Q simply 
by replacying G and G' with G and G' . These formulas are well-suited to computer 
computations, so in later sections we shall use G and G' as needed on graphs with one more 
cut edges without showing the steps used to compute these quantities. 

8 Loop-erased random walk 

In this section we compute the probability that the LERW in Z^ from (0, 0) to oo passes 
through the vertex (1,0), and the analogous probabilities for other points and for other 
lattices. 

8.1 General remarks 

We let Py^w denote the probability that the LERW started from (0, 0) to oo passes through 
edge {v,w), in the direction from v to w . Likewise we let P^ denote the probability that 
the LERW passes through vertex w . It is straightforward that 

Our strategy is compute these edge probabilities. 

We find it conceptually convenient to work with finite graphs Q, set up our equations 
for spanning tree and grove event probabilities in terms of the finite-graph Green's function 
G^ and its derivative (G^)' using the formulas from section [6], and then afterwards take the 
limit ^ — !■ L using the formulas from section [71 We are at liberty to use any convenient 
sequence Qk of finite graphs that converge to the lattice L, since the limiting measure on 
spanning trees of L is independent of the choice of Qk , and the event that the LERW from 
(0, 0) to oo uses a given edge is a measureable event in the limiting measure. So we choose 
a sequence Qk for which it is convenient to compute {G^'')' , as described in section § 7.2 



These graphs Qk have a wired boundary vertex (see [Figure 3 ) that we will label oo, even 



though the graphs are finite. The other vertices of Qk we will label by their coordinates in 
L. We define P^^ and P^ in the same way as we defined P^^w and P^; for fixed v and w, 
limg_>.L P„ ,„, = Pv,w ■ 

We remark that once an edge traversal probability P^^ has been computed, finding the 
corresponding probability P^,^ for the reversed edge is straightforward: 

Lemma 8.1 ( [KenOObj ). 
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Proof. Let the origin be at the center of a large box with wired boundary whose size we 
will send to infinity. For simple random walk started at , the expected number of traversals 
of {v,w) minus the expected number of traversals of {w,v) is the edge conductance Cv,w 
times the difference in Green's functions. The same also holds for loop-erased random walk 
started at 0, since cycles are reversible. D 



The intensity of the undirected edge {v,w} is -Pj^^i 



pL 



P^^, and the undirected 



edge intensities turn out to be nicer numbers than the directed edge intensities. The vertex 
intensities are easily calculated from the undirected edge intensities, and given the potential 
kernel of L, the directed edge intensities are easily recovered from the undirected intensities. 
We show below how to calculate the specific edge intensities relevant to the -P(i,o) com- 
putation before describing the general method that can be used for any edge. 



8.2 Loop-erased random walk on Z 

We first recall the potential kernel for Z^ in Figure 10 
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Figure 10: Square lattice potential kernel. 
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8.2.1 Probability LERW uses edge (1, 1)(1, 0) 



To compute P(i,i),(i,o) , the probability that the LERW uses 
edge (1, 1)(1,0), in fact we use the annular surface graph shown 
in Figure 3, and shown more schematically on the right. There 
are four nodes, whose coordinates are (1,0), (1,1), (0,0), and 
oo, and which we also refer to as nodes 1, 2, 3, and 4. 

Consider a uniformly random spanning tree of this graph Q , 
and within it the path from (0,0) to oo, i.e., the path from 
node 3 to node 4. Those spanning trees in which the path from 
(0, 0) to oo contains the directed edge (1, 1)(1, 0) are in bijective 
correspondence with groves with connection type 3, 2 1 1,4; the bijection is simply to delete 
the edge (1, 1)(1,0) from the spanning tree, or in the reverse direction, to add the edge to 
the grove. Therefore, 

p5 — ^[3,2|1,4] _ ^^,g ^ ^,g ^ ^,g 




(i,i),(i,o) 



r"g ri'Q ri'Q i r^Q r'Q 

-Lr ^2 - ^ 2,3~ ^ 3,1 + *-^2,3 ~ *-^l,J 



Z[l,2,3,4] 

by (llOcj) from section O Using the method described in section [71 we (or the computer) can 
compute 





G^' 


(1,0) (1,1) (0,0) 


G'^' 


(1,0) (1,1) (0,0) 






(1,0) 


-i -1 

"J 4 4 


(1,0) 


-^ -^ 

" 32 32 






(1,1) 


_i -i 

4 ^ 7r 


(1,1) 


32 ^ 27r 






(0,0) 


-i -- 

4 TT 


(0,0) 


.32 2tt ^ . 






We can now evaluate Vmvg^ip. PH -^^ ^ g^ by substituting G^ for G^ = G^ and G 


G' 


s = {Gsy-. 

^2 3 15 11 3 1 

^(i,i),(i,o) - 32 + 2^ 32 ^ + 4 -+16 27r 






pZ2 _ 1 
^{(1,0),(1,1)} - o- 









for 



(23) 



With P^ ^s ^ QN in hand, it is possible to compute P^ qn using an ad hoc argument and 
existing techniques |BP93| IKen97] for computing finite-dimensional marginals of the random 
spanning tree. However, here we take the systematic approach of computing the probabilities 
that various edges are on the LERW path. 



8.2.2 Probability LERW uses edge (2,0)(1,0) 

Next we show here how to compute P^ q. ^ qn using the methods developed in sections [6] 
and [3 The vertices of interest are (0,0), (2,0), (1,0), and oo, so we make them the nodes 
of an annular-one graph. Since vertices (0, 0) and (2, 0) are not on the boundary of the same 
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face in Z^, we cut (undirected) edge (1,0) (1, 1) of the lattice, so that all the nodes (other 
than oo) are on the boundary of the same face. The endpoints of any cut edge also become 
vertices of interest, so they are included as nodes in the resulting annular graph. Thus our 
annular graph has five nodes: (1, 0) , (2, 0) , (1, 1) , (0, 0) , and cxd, which we also refer to as 
nodes 1, 2, 3, 4, and 5. Nodes 1,2,3,4 are arranged counterclockwise on the inner face 
with a zipper between nodes 1 and 4, as in section [61 

Next we relate the spanning trees of the original (uncut) graph to groves in the graph 
without (undirected) edge {(1, 0), (1, 1)}. Suppose the LERW from (0,0) to oo passes 
through edge (2, 0)(1, 0) (which we also call edge 2 — )■ 1). There are three cases depending 
on whether or not edge (1, 1)(1,0) or (1,0)(1, 1) (edge 3 — )> 1 or 1 — )■ 3) also occurs in the 
tree: 




1. Edge 3 — 7- 1 occurs in the tree. Trees in this case are in bijective correspondence with 
groves of the cut graph of type 4, 2|3|1, 5, where the bijection deletes edges 2 — )■ 1 and 
3—7-1 from the tree to get the grove, or adds them to the grove to get the tree. 

2. Edge 1 — 7- 3 occurs in the tree. Trees in this case are in bijective correspondence with 
groves of the cut graph of type 4, 2|1|3, 5 (with the obvious bijection). 

3. Neither 3 — )■ 1 nor 1 — ?► 3 occurs in the tree. Trees in this case are in bijective 
correspondence with groves of type 4, 2 1 1,5 (with the obvious bijection). (When a 
node is not listed in the partition, it is treated as internal, and could in principle occur 
in any of the parts, though in this case it would necessarily occur in the 4, 2 part.) 



Putting these three cases together, we obtain 

Z^[4,2|3|l,5] + Z^[4,2|l|3,5] + Z^[4,2|l,5] 



pQ 

M2,0),(l,0) 



pG 

M2,0),(l,0) 



Z^[4,2|3|l,5] , Z^[4,2|l|3,5] Z^[4,2|l,5] 



Z^ 



ZG 



ZG 



X 



Z^ 
ZG 
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Each of the expressions within the brackets is of a type which section |6] shows us how to 
compute in terms of G^ and G . These are 



Z[4,2|3|l,5] 



z 

Z[4,2 13 


5] 


Z 

Z[4,2 1 


5 



1,3^^2,3 + tjl,3*-^2,3 — tjl,4tJ2,3 + tjl,4L^2,3 + '-^l,3'-^2,4 ~ '-^l,3'-^2,4 

l,2'-^3,3 + '-^l,4'-^3,3 ~ '-^l,4<-^3,3 ~ '-^2,4'-^3,3 + '-J2,4'-J3,3 ~ '-Jl,3'-J3 

+(^2, 36*3 4 + G*! 2^3,4 — G*]^ 36*3,4 + G'l,3G'3,4 + G'2, 3^3,4 — 6*2, 3^*3,4 

1,2'-^1,3 ~ '-^1,2'-^1,3 " '-Jl,2'-^l,4 + '-^1,3'-^1,4 + '-^1,2'-^1,4 ~ '-^1,3'-^1,4 
l,l'-^2,3 + '-^l,4'-^2,3 + '-^l,l'-^2,4 ^ '-^l,3'-^2,4 ~ ^\,\^ZA + '-^l,2<-^3,4 

-G-^ 2 + ^1 4 ~ ^1,4 ~ G2 4 + Cr2,4 



From section [7] we know how to compute G^ and G"^ and also the ratio Z^ /Z^ in the hmit 
Q ^ 7? . These values are 



lim — -T 

g^i? Z^ 



1 

2' 



and 



^Z2\{{1,0)(1,1)} 

(1,0) 
(2,0) 

(1,1) 
(0,0) 



and 



(5«2\{(1,0)(1,1)} 

(1,0) 

(2,0) 

(1,1) 
(0,0) 



(1,0) 
1 

8 
1 3 

2-K 8 
3 



27r 



(1,0) 



J ]_ 

32 167r 

_3_ 

16 

_9 T_ 

32 167r 



(2,0) 



J 

27r 
+ 77 



27r 



.7 + A+ 1 

S T^ 7r2 T^ ,r 



(2,0) 



(1,1) 

_3 

8 

1 3_ 

8 27r 

1 

8 

1 3_ 

8 27r 



;i,i) 



1 

167r 



_3_ 
32 



_3_ 

16 



16 



_9 1_ 

167r 32 

_ 3 _ 

2^ 



_1_ 
32 



9 

167r 



9 

167r 



J_ 
32 



Upon substituting and simplifying we obtain 
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(0,0) 



1 

27r 


3 

8 


+ 


2 , 1 


1 

8 ' 


3 

27r 


-f 


2 1 




7r2 TT J 



(0,0) 



7 
167r 



_9_ 
32 



9 I 3 I 5 
16 ~'~ 27r2 "■" Stt 



32 



9 

167r 



M(i 



{(1,0),(2,0)} 



(24) 
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8.2.3 Probability LERW passes through (1,0) (sandpile value) 

Combining ( J23l) and (!24l) together with the symmetry relations P^ _^-. ,^^ q-. 



P, 



(0,0), (1,0) 



1/4, we obtain 



Mi,o) - T^ 



16 



-^(1,1), (1,0) ^^"^ 



(25) 



8.2.4 Probability LERW uses edge (2, 1)(1, 1) 

We delete edge (1,0) (1, 1) from Z^'s finite graph approximation Q , the same edge that we 
deleted in section |H21 and again view Q = Q \ (1,0) (1, 1) as an annulus with 5 nodes, but 
this time nodes 1-5 are (1,0), (2,1), (1,1), (0,0), and oo. 

For spanning trees that use edge (2,1) (1,1), there are again three cases depending on 
whether or not and in what direction the cut edge {(1, 0), (1, 1)} is used. Spanning trees 
in each of these cases are as before in bijective correspondence with groves pairings on the 
annulus: 




At this point it should be clear how to carry out the computation, which we did by 
computer, so we just state the answer: 



^(2,l)(l,l) - +T7 + 



16 47r2 



in 



M(i,i),(2,i 



1 



1 



8 27r2 Air' 



(26) 



8.2.5 Probability LERW passes through (1, 1) 



Since vertex (1, 1) lies on a symmetry axis, two edge probabilities suffice to determine P^ ^^ . 
Combining (!26|) and (!23|) we obtain 



^(1,1) 



- + At- — = 0.221083... 
4 27r2 An 



(27) 



This compares well with the empirical probability of 0.221088 based on 10^ LERW's from 
the center of a 1400 x 1400 grid with wired boundary. 
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8.2.6 Probability LERW uses edge (2, 1)(2, 0) 

Again we delete edge (1, 0)(1, 1) from l? 's finite graph approximation Q . This time we need 
six nodes, which in order are (1,0), (2,0), (2,1), (1,1), (0,0), and oo. There are three 
main cases as before, depending on whether or not and in what direction the spanning tree 
of Q uses edge (1,0)(1, 1) (which we now refer to as 1 — )■ 4 or 4 — )■ 1), but now there are 
several subcases. The edge we are studying is 3 — >■ 2 : 




Once we have this list of cut edges, nodes, and grove types, we can carry out the computation 
and obtain 



M2,l 



1 



(2,1)(2,0) 



^TT^ 47r 



{(2,0),(2,1)} 



47r2 27r 



(28) 



8.3 General edge intensities 



The calculations in section 18.21 illustrate the general method for computing the probability 
that LERW passes through an arbitrary given edge. We identify a set of edges to cut so as 
to place the starting point of the LERW and the endpoints of the edge on the same face. 
These three vertices, together with the endpoints of all of the cut edges, comprise the nodes 
on the inner boundary of the annulus. We number these nodes in counterclockwise order so 
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Figure 11: Undirected edge intensities of loop-erased random walk on Z^. The edge- 
intensities of {x,x){x,x — 1) and {x,x — l){x + l,x — 1) are identical for a; = 1,2,3 (and 
perhaps larger values), despite there being no lattice symmetry that would imply this. 
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that the zipper starts between nodes 1 and n — 1. The vertex labeled oo becomes node n, 
which is on the other boundary of the annulus. For each possible set of the cut edges that 
might occur within a spanning tree of the original graph Q , we associate spanning trees 
containing precisely those cut edges with groves in the cut graph Q. We can enumerate 
all possible compatible grove partitions, and then rewrite these as a linear combination of 
partial pairings as described in section [6^ Using the methods described in sections [6] and [71 
we compute the probability that the resulting grove is of the appropriate type for the edge 
of interest to be on the path from (0, 0) to cxd. 



8.4 Loop-erased random walk on the triangular lattice 

In this section we sketch the computation showing that on the triangular lattice, a LERW 
from (0, 0) to oo passes through vertex (1, 0) with probability 5/18. (We choose coordinates 
for which vertex {x,y) is the one located at position {x + y/2,y\/3/2) in the plane.) The 



starting point of the calculations is of course the potential kernel, which is shown in Figure 12 



The undirected edge intensities are shown in Figure 13 
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Figure 12: Triangular lattice potential kernel. 

The following diagrams indicate how to compute the probability that, respectively, edge 
(0,1)(1,0), edge (1,1)(1,0) and edge (2,0)(1,0) are on the LERW from (0,0) to 00. The 
methods of sections [7] and [HI automate these computations. 
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ptriangular 

^(0,1),(1,0) 
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{(1,0),(0,1)} 



1 
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ptriangular 

Mi,i),(i,o) 
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triangular 



{(1,0),(1,1)} 



31 

108 
2 

27' 



73 

27r 




^triangular V <J '^•^ 

(2,o),(i,o) - V " 108 



itriangular 



Combining these probabilities we obtain the "sandpile value" of p'^^^^^^ ^^ 



{(1,0), (2,0)} 27' 

^trian 
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Figure 13: Undirected edge intensity of loop-erased random walk on the triangular lattice. 
The edge-intensities of (x, x){x, x — 1) and {x,x — l){x + l,x — 1) are identical for a; = 1, 2 
(and perhaps larger values), despite there being no lattice symmetry that would imply this. 
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8.5 Loop-erased random walk on the hexagonal lattice 



The potential kernel of the hexagonal lattice is shown in [Figure 14) Using the methods 
described above, we computed the undirected edge intensities for a LERW from (0, 0) to oo, 



which are shown in Figure 15 
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Figure 14: Hexagonal lattice potential kernel. 



8.6 Loop-erased random walk on Z x M 

We consider next a weighted version of 7? , where each horizontal edge has weight c, and each 
vertical edge has weight 1/c. This graph is isoradial, so we may compute the Green's function 
using |Ken02] . Because the lattice is symmetric under a 180° rotation and invariant under 
translations, we can also compute G' . This gives us all the necessary information we need 
to compute the probability that LERW from (0, 0) passes through (1, 0) . It is convenient to 
let c = tan^. After a computation similar to the ones above, we find that the LERW passes 
through vertex (1,0) with probability 



1 9 

4^2^ 



02 



TT^ sin^ I 



?)■ 



When 9 = tt/A, we have c = 1, and this above probability reduces to 5/16, in agreement 
with our earlier calculation for Z^ . 



53 































91 13 22 


/ V 




'A* 
aV 


A'" 

„_, ... ^.. ... '/; 






MS 6.' ,/3-„ 


%\^ 




/:4' 


\ • > 

A*- 

a; * 




/ 


5 


"4" 


167 11 17 


•A 

A 

A 




\ 






36 












< 




-X 


5 


% 


A* 


.'4" 


36 




„^.,...,...,^ 


o.^^ 




X- 




^>* 


. :4' 


'~4« 


(origin) 


1 








5 7 -V,, 


4^ 


3 




81 36^fJn 





Figure 15: Undirected edge intensity of loop-erased random walk on the hexagonal lattice. 
Using coordinates for which vertex {x,y) is the one located at position (x — y/2,y\/3/2), 
the edge intensities (3x — 1, 3x — l)(3x — 1, 3a; — 2) and (3a; — 1, 3a; — 2)(3x, 3a; — 2) are 
identical, as are the intensities for edges (3a;, l)(3a; + 1,2) and (3a;, l)(3a;, 0) (for x = 1,2 
and perhaps larger values), despite there being no lattice symmetry that would imply this. 
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Figure 16: Portion of uniform spanning tree and LERW from (0,0) to oo on Z x (Z/20), 
which approximates Z x M. The probabihty that the LERW from (0,0) passes through 
(1,0) is l/4-l/7r2. 



In the isoradial embedding of the lattice into the plane, if the horizonal edges have 
length 1, then the vertical edges have length c. An interesting special case is the limit 
c — )■ 0. Then random walk on this weighted graph converges to a standard Brownian 
motion in the vertical direction, except at a Poisson set of times with intensity 1 , where the 
walk jumps left or right with equal probability. The random walk on this graph is then a 
continuous-time random walk on Z in the horizontal direction and a Brownian motion on 



M in the vertical direction. Figure 16 illustrates a portion of the uniform spanning tree on 
this graph, together with the path from (0,0) to oo. From the above formula, we see that 
in this limit, the probability that the LERW passes through (1, 0) converges to 1/4 — I/tt^. 



A The annular-one matrix A 



n 



A.l Association of rows and columns 

Recall the matrix A^ for the annulus with n — 1 nodes on one boundary and 1 node on the 
other boundary. We prove that A„ is invertible. More specifically, we prove that det A^ is 
a polynomial of degree 2"~^ — \{J^i^ with lead coefficient ±1 and constant coefficient ±1. 
A key part of the proof is a natural correspondence between the rows and columns of 
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the matrix A„. This correspondence is based on the "cycle lemma" bijection of Dvoretzky 
and Motzkin |DM47] (see also |DZ90] ). This bijection maps a cyclic sequences oi k + 1 +'s 
and k — 's to a Dyck path of length 2k together with an "offset", a number between 1 and 
2/c + 1. (There is a unique + that can be deleted so that the subsequent sequence defines a 
Dyck path; the position of this + is the "offset". Given the Dyck path and offset, the Dyck 
path is interpreted as a sequence of +'s and — 's, a + is prepended to the this sequence, and 
the offset determines a cyclic shift.) In our setting, k = n/2 — 1 , in det J^^ , the indices in R 
give the locations of the +'s, and the indices in S give the locations of the — 's. The offset 
that the cycle-lemma bijection gives (a number from 1 to n — 1) determines the index with 
which n is paired, and the chords under the Dyck path determine the rest of the pairing. 
The following example illustrates the bijection adapted to (n — 1, 1) -annular pairings: 



,8^4 



^8^9, 



/7'^ I'O h' ,8^, 



{1,2,5,9,10,12} / khA 4^y I y KyA 5 2|10|9|1|12 

{3,4,6,7,8,11} ^ 'v / V .Vr! ^ 4-^^e ^^ 4l6l 7 Isllll 3 

f^. ^ 1"2 

and in the reverse direction, 



1,1112, 613, 1214,517,1018,9 




5|2|10|9| 1 |12 , 5,2,10,9,1, 12 



, xx|^, w|o, x^|-*, ^1 I , xw|u, J ^ 3 4 5 6 7 8 9 I'O 11 1 2 3 ^ 4l6l 7 ISlllI 3 ^ 4,6, 7, 8,11, 3" 



A. 2 Invert ibility 

We prove more than lTheorem 6.Ti 

Theorem A.l. The determinant of the annular-one matrix, det An, as a polynomial in (, 
takes the value 1 when C = 0, has degree 2"^^ — K^'^g) ' ^^'^ ^^^ leading coefficient ±1. 

Proof. It is convenient to work with a slightly different annular matrix A* which corresponds 
to a symmetrized placement of the zipper, with a zipper starting from each interval of the 
first boundary between each pair of nodes. We do this by splitting the original zipper into 
n — \ zippers each with parallel transport 2;^/("~^) , and then deforming these zippers so 
that their endpoints lie in each of the n — 1 intervals between the nodes. When we deform 
a zipper across node i{^ n) in the counterclockwise direction, then the parallel transport 
from i to any other node j is multiplied by ^^Z^""^) . For each column of the annular matrix, 
say indexed by directed pairing a, the column is scaled by z^^'^"^ '' according to whether 
node i is a source or destination in a . Likewise, each row of the annular matrix, say indexed 
by det=SfJ^, is scaled by 2;^i/("~i) according to whether i E R oi i E S . Thus the effect 
of deforming these zippers is to conjugate the annular matrix by a diagonal matrix, so in 
particular det A* = det A„ . We change variables to 

^ ^ ^S/Cn-l) ^ ^l/(n-l) 

56 



so that the nonzero entries of A* are integral powers of w. For example, the first two rows 



of Ag are 



^lO 


coin 


to-!* 


CO^ 


;OCO 


;OCO 


^CS 


CO(M 


CO^ 


to^ 


'^CO 


^^ 


CO(M 


coin 


CM^ 


IM^ 


^in 


^co 


LO^ 


LOC<l 


IM^ 


COIM 


^lO 


(M^ 


LO^ 


^lO 


^co 


lO-^ 


COIM 


-5t<CO 





i^ 


;^ 


;^ 


;^ 


i^ 


i^ 


;^ 


;^ 


'M 


;^ 


det^?:i^ 


1 


—w 








—w 











w^ 


—w 


det^lt;^ 





1 

















w^ 


—w 






and the other rows are determined by cyclic rotations. 

From the determinant formula, we see that the diagonal entries of A* are all 1 . Consider 
the column indexed by directed pairing a. Since A* is symmetric under cyclic rotations of 
the indices 1 , . . . , n — 1 , let us assume for convenience that a pairs n — 1 to n , so that we 
can write 



a 



0,1,1 I 
ai,0 I 



lO-n/2-1,1 I n 

\o-n/2-i,a In— 1' 



Referring to the above bijection, since n — 1 pairs to n, for each j we have Oj^q < 0'j,i- 
Column a contains 2"/^"^ nonzero entries, one for each sequence /i, . . . , fn/2-1 of n/2 — 1 



/2 

O's and I's, where the row is indexed by 

d t ^"'^'^'^' ' " ' ' ""/^-^•^"■^"/a-i ' '^ 
ai,/i,---,an/2-i,/„/2-i''^ ~" 1 ' 

Since the pair {aj^, aj^i) crosses a^^i — aj^ zippers, this pair contributes W'^^ to the matrix 
entry. In particular, all nonzero nondiagonal entries of A* have positive powers of w. This 
implies 



det A^ 



detA! 



which of course implies that the annular matrix is nonsingular. 

We can also compute the degree of the polynomial det A* . For a given column, the row 
that maximizes the power of w is the one for which fo, ■ ■ ■ , fn/2-1 = 1, • • • , 1, and the power 
is the area under the Dyck path. The mapping from a column a to the row det =Sf^ which 
has the highest power of w is also a bijection, in fact it is a simple variant of the bijection 
described above. This implies that the leading coefficient of the polynomial detA* is ±1, 
and the degree is 

deg det A* = (n — 1) x N^ area under Dyck path. 

Dyck paths of length n — 2 

For Dyck paths of length 2k = n — 2, the above sum is known (see Sloane's A008549) to be 



4^= _ 



2k + l\^^^_,_fn-l \ ^^n-2_\( n 
k ) \n/2 -l) 2 \n/2 

Consequently, the degree of the polynomial det A„ (a polynomial in ^ ) is 2"~^ 



-lU)- n 
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A. 3 Cancellation of higher order derivatives 

Next we restate and prove ITheorem 6.31 

Theorem A. 2. For all positive even n and pairings o o/ {1, . . . , n}, Z„jZx\2\-\n is a poly- 
nomial of degree n/2 in the quantities 

{Lij ■.l<i<j<n} and {L'-j : I < i < j < n - I}. 

Proof of \Theorem 6.3\ . For a general finite graph with general parallel transports, define 
^[cr] as in the case of a planar graph or annular-one graph. Let us then define D^ to be 
the result of multiplying row {R, S) of the matrix A„ by the vector of ^[aj's. For annular 
graphs we have D^ = det=SfJ^, but for general graphs these two quantities will be different. 
We can "recover" the ^[cr]'s from these -D^'s by multiplying by A~^: 

^[1|2|^--H = (1 - Qm. E «<^-^.^(C)^i (29) 

-TtjO 

where 

and m^ > is the smallest integer such that for each {R, S) , the coefficient {l—Cy'^'^{/^n^)a,R,s 
is nonsingular at C = 1 • 

Let us consider now the complete graph on n nodes, so that ^ = —A and J^ij is the 
edge weight between nodes i and j , times the parallel transport to i from j . Note that in 
this case =Sf is (except for the diagonal entries) a general Hermitian matrix. Each Df^ is a 
polynomial in the entries of ^ , with coefficients that involve powers of ( . (For the complete 
graph, any matrix times the vector of ^[cr]'s will yield polynomials in .^.) 

Next we change variables by setting ( = e^ . This makes it more convenient to differentiate 
the annulus response matrix, since ^^j,i = ~'^'^i,j- The C ~^" 1 limit is of course equivalent 
to t — !■ 0, and (1 — e*)*" has a zero of order m at t = 0, with ^(1 — e*)™ = (—l)"^m\. 

For general nonzero edge weights and smooth (in t) parallel transports on the complete 
graph, ^[a] is finite and 5°[1|2| ■ ■ ■ |ra] = 1, so for any i < m„ it must be that we get zero 
when we differentiate the numerator from fl29l) . i.e., YliRs^'^Asi.OD^, i times with respect 
to t and then set t to 0: 



o-E E 



d'-'^ 



k=0 ^ ' R,S 



O^a^R.S 



d 



k 



t=0 



dt^^"^ 



t=0 

k 



Since J^ is generic, we can rescale the kth derivative of each J^ij by a factor of (3 , and 
deduce that for each k, i with k < i < m„ 

o = E;£^"-«.-,_„x:l^«H,_„- (30) 



«,. *'- 



t=0 
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Next we write -D^() for the polynomial function of a Hermitian matrix, which, when 
evaluated on the response matrix =Sf of the complete graph, gives Df^ = D^{^) . Let 
di denote the differential operator for which diJ^ij = ^-^-^ij and di^j^i = \-^-^j,i but 
diJ^h,j = for h,j 7^ i. For the complete graph, each monomial of the polynomial D^^ 
includes each index i exactly once. (This also holds for annular graphs, though of course 
the polynomials are different.) Using this property of these polynomials, we can write the 
kth derivative of ^ = =^(t) as follows: 



Jk 



J2 dt,---d,^Dl,{^). 



«i,«2,---G{1, .-■,"} 



(31) 



Given a set T oi k nodes, for each i E T and each j we rescale 4=Sf^, ,• by a factor 

at 'J ^^Q 

of /3, without changing any of the other derivatives at t = 0. If i, j G T, then we rescale 

■jf2-^ij by a factor of /3^. (Recall that ^ is generic Hermitian, so we can do this.) 
dt 'J ^^Q 

The coefficient of P^ within the fcth derivative is obtained from f l3T|) by including only those 

terms for which ii, . . . ,1^ is a permutation of T. Then substituting fl3T]) into fl30|) and taking 

the coefficient of /3^, we find 



= 5Z 'lIF^.^'^'R^s 



R,S 



dt^- 



t=0 



D^(^(^)) 



i=0 



(32) 



where J^^'^^ where is the Hermitian matrix obtained from =Sf by replacing Jtfij with ^^i,j 
for each i E T and j ^ T oi i ^ T and j E T, and replacing J^ij with -^^J^ij for each 

From our definition of D^ , for the complete graph we have 

where the sum is over directed pairings r that happen to be annular. Next we take the 
pairing a , and for each pair {i,j} of a, we rescale =Sfjj by a factor /3. Then the coefficient 
of /S""/^ in ( 132|) only arises when r = a in the above sum, so 



E 

R,S 



dt^- 



:Ota,R,S 



t=Q 



X {An)R,S,r X n "^* 



(T) 



ieRJeS 



t-0 



(33) 



The product term is the formula always takes the same (generically nonzero) value, except 
for a sign, which is given by the parity of the number of indices in T which are also in S . 
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So we cancel this factor (keeping the sign), rewrite a in terms of A„, and obtain 







R,S.cT pairs R to S 



[iA-\e')UnA^ - eT^] 



_l\\SnT\ 



t=0 



(34) 



whenever \T\ < i < m. This is purely an identity about the matrix A„^. 

Next we return to annular graphs, and differentiate the numerator m = ma times: 



fc=0 






R,S 



dt"-- 



X -- det ^^ 



R 



t=0 



We can rewrite the kth derivative of det=SfJ^ using the differential operators di, as in ( l3T]l . 
Let T be the set of nodes for which we applied di an odd number of times. If for some i 
we applied di more than once, then the size of the set T will be less than k, and then from 
( 134|) we see that the coefficient of such terms is . 

Next, suppose the variable J^i^n is differentiated. Since n is in each set 5", we may as 
well replace T with T\ {n} and introduce a global sign, but then since \T\ is smaller, we 
see from (1341) that the coefficient of such terms is . D 
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